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A  mathematical  model  of  the  closed  solar  air  heaters  is  used  in  particular,  to  assist  in  interpreting  the 
observed  phenomena  in  the  solar  air  heaters,  to  design  the  system,  to  predict  the  trends,  and  to  assist  in 
optimization.  In  this  paper,  various  mathematical  models,  mainly  analyzing  the  heat  transfer  process  of 
solar  air  heaters,  are  reviewed  and  classified  based  on  the  model,  the  number  of  the  cover,  the  shape  of 
the  absorber  and  the  presence  or  not  of  the  packing  bed.  Although  the  models  have  evolved  to  a  point 
where  several  features  of  the  process  can  be  predicted,  more  effort  is  required  before  the  models  can  be 
applied  to  define  actual  operating  conditions  as  well  as  to  further  investigate  new  closed  solar  air 
heaters.  It  is  shown  that  the  major  governing  equations  in  the  models  are  based  on  the  first  law  of 
thermodynamics. 
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1.  Introduction 

One  of  the  most  potential  applications  of  solar  energy  is  the 
supply  of  hot  air  for  the  drying  of  agricultural,  textile  and  marine 
products,  and  heating  of  buildings  to  maintain  a  comfortable 
environment  especially  in  the  winter  season.  Several  designs  for 


Nomenclature 

A 

area  element 

B 

height  of  solar  air  collector 

C 

specific  heat 

D 

depth  of  air  channel 

Dh 

equivalent  diameter  of  air  channel 

h 

heat  transfer  coefficient 

k 

thermal  conductivity 

L 

Length  of  air  heater 

M 

mass  per  unit  collector  area 

m 

mass  flow  rate 

Nu 

Nusselt  number 

Pr 

Prandtl  number 

R 

recycle  ratio 

r 

fraction  of  the  mass  flow  rate 

Ra 

Rayleigh  number 

Re 

Reynolds  number 

S 

incident  solar  radiation 

T 

temperature 

t 

time 

U 

loss  coefficient 

V 

wind  speed 

W 

width  of  absorber  (collector  width) 

x,y 

space  coordinate 

Greek  letters 

a 

solar  absorbance 

S 

thickness 

£ 

emissivity 

a 

Stefan-Boltzmann  constant 

X 

Transmitivity 

P 

reflectance 

Pf 

fluid  density 

Subscript 

a 

ambient 

b 

bottom  plate 

c 

convection 

/ 

fluid  (air) 

g 

cover 

i 

inlet 

0 

outlet 

P 

absorber  plate 

r 

radiatif 

T,t 

top 

solar  air  heaters  have  been  proposed  and  discussed  in  literature. 
The  designer  and  potential  user  of  these  systems  must  consider  a 
number  of  factors  when  comparing  their  merits.  These  can  mainly 
be  categorized  as:  (i)  thermal  performance,  (ii)  cost  and  (iii) 
lifetime,  durability,  maintenance  and  ease  of  installation.  Thermal 
performance  of  collectors  is  compared  by  using  the  concept  of 
thermal  efficiency.  It  is  generally  believed  that  the  thermal 
efficiency  of  a  solar  collector  is  the  major  requirement  for  the 
prediction  of  thermal  performance  of  the  complete  solar  system  of 
which  the  solar  air  collector  is  a  part  [1  ].  It  is,  therefore,  important 
that  this  thermal  efficiency  information  be  in  a  form  which  is 
useful  for  existing  and  projected  future  solar  energy  system  design 
methods.  To  save  energy,  the  method  of  exergy  analysis  is  useful  in 
solar  systems,  because  it  appears  to  be  an  essential  tool  for  system 
design,  analysis  and  optimization  of  thermal  systems  [2].  One 
indication  is  that  there  are  about  two  review  papers  in  the  area  of 
solar  air  heaters  published  in  recent  years.  The  first  was  proposed 
by  Chandra  and  Sodha  [1].  These  authors  have  provided  a 
fundamental  understanding  of  testing  procedures  for  solar  air 
heaters.  The  second  was  proposed  by  Ekechukwu  and  Norton  [3]. 
They  have  classified  solar  air  heaters  broadly  into  two  types:  bare 
plate  and  cover  plate  solar  energy  collectors.  Based  on  this 
classification,  these  authors  have  summarized  various  designs  of 
solar  air  heaters. 

With  the  development  of  computer,  hardware  and  numerical 
methodology,  advanced  mathematical  models  are  being  used  to 
carry  out  critical  investigations  on  solar  air  heaters.  The  advantages 
of  these  methods  are  that  they  can  produce  extremely  large 
volumes  of  results  at  virtually  no  added  expense  and  it  is  very 
cheap  to  perform  parametric  studies  to  optimize  equipment 
performance.  The  second  reason  for  such  work  on  numerical 
simulation  is  that  some  parameters  are  difficult  to  test,  and 
experimental  study  is  expensive  as  well  as  time  consuming. 

The  purpose  of  this  work  is  to  review  the  present  state  of 
mathematical  modeling  of  solar  air  heaters.  The  validation  is  an 
important  step  in  mathematical  modeling  development,  and 
therefore  comparisons  with  actual  experimental  values  or  theore¬ 
tical  results  have  been  included  where  possible.  Since  most  of  the 
models  have  to  be  solved  numerically,  the  numerical  techniques 
employed  in  the  solution  of  different  models  are  presented. 

2.  Mathematical  models 

2.2.  A  bare  plate  solar  air  heater 

This  consists  of  a  single  channel  design  with  single  air  flow 
between  absorber  and  bottom  plates  with  insulation  (see  Fig.  1). 
Choudhury  et  al.  [4]  studied  the  thermal  performance  of  this 
model.  The  energy  balance  equations  for  the  absorber  plate  and  the 
rear  plate  have  been  written  as 


Bare  plate  :  apS  =  hpa(J P-Ta)  +  hpb(T P-Tb)  +  hpf(T P  —  Tf) 

(1) 

Rear  plate  :  hpb(Tp  -  Tb)  =  hbf(Jb  —  Tf)  +  Ub(Tb  -  Ta)  (2) 

With  a  design  selection  methodology  whereby,  for  combined  fixed 
values  of  pressure  drop,  air-flow  rate  and  duct  length,  duct-depth 
values  and  the  resultant  values  of  forced  convective  heat-transfer 
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Fig.  1.  A  schematic  view  of  a  bare  plate  solar  air  heater  [4,7-9]. 
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Fig.  2.  A  schematic  view  of  a  single  cover  solar  air  heater  [7,9,13,16]. 


coefficients,  air-temperature  increment  per  unit  incident  flux  and 
system  efficiency  have  been  computed.  From  the  discussion,  it  was 
shown  that  the  computation  methodology  and  the  design  curves 
obtained  can  be  successfully  used  to  construct  solar  air  heaters 
with  pre-determined  acceptable  values  of  pressure  drop  with  the 
required  rate  of  air  flow,  minimum  material  cost  and  highest 
possible  efficiency. 

Thermal  performance  of  this  design  has  also  been  predicted  by 
Ong  [7].  The  mathematical  model  proposed  by  this  author  is 
similar  to  the  one  of  Choudhury  et  al.  [4]  with  a  small  change  in 
Eq.  (1 )  where,  hpa  is  replaced  by  Ub  =  hw  +  hrgs.  In  this  mathematical 
model,  the  solar  collector  was  assumed  sufficiently  short  for  which 
the  assumptions  were  valid.  The  mathematical  solution  procedure 
involved  a  matrix  inversion  of  the  mean  temperature  vector 
derived  from  energy  equations.  Predicted  temperatures  for  the 
solar  air  heaters  were  presented.  The  effects  of  wind  and  film  heat 
transfer  coefficients  on  the  prediction  were  discussed. 

Njomo  [8]  and  Njomo  and  Daguenet  [9]  have  investigated  heat 
transfer  in  the  above  design  (Fig.  1).  With  general  assumptions  in 
the  heat  transfer  modeling,  the  energy  balance  equations  of 
different  components  of  the  solar  air  heater  have  been  written  as: 

dT 

For  selective  absorber  :  MPCPP  =  apS  -  hrpb(Tp  -Tb) 

—  h-cpf(Tp  —  T f)  —  (hcpa  +  hr pa) (T p  —  T a)  (3) 


c  •  fi  n  ^Tp  rrifCpf  dT p 

For  air  flow:  pfefCpf- + 

=  hcpf{Tp  -Tf)  +  hcbf(Tb  -Tf)  (4) 

dT 

For  bottom  plate  :  MbCpb 

=  hrpb(Tp  -  Tb)  -  hcbf{Tb  -Tf)-  Ub(Tb  -  To)  (5) 

With  the  mathematical  model  proposed,  Njomo  [8]  predicted 
the  thermal  performance  of  such  a  collector.  To  determine  a  solar 
air  heater  thermal  behavior,  Njomo  and  Daguenet  [9]  used  a 
sensitivity  analysis  which  is  an  excellent  technique  to  help 
designers  and  decision  makers  in  anticipating  and  preparing  for 
a  solar  system  project.  The  results  have  demonstrated  the 
existence  of  an  important  dimensionless  parameter  referred  to 
as  thermal  performance  factor  of  the  collector  that  compares  the 
useful  energy  that  can  be  extracted  from  the  heater  to  its  overall 
heat  losses.  A  detailed  sensitivity  analysis  has  been  performed,  and 
the  results  successfully  summarized  to  measure  the  impact  of 
changes  in  operation  and  meteorological  parameters  on  the 
simulation  values  of  the  thermal  efficiency  and  the  fluid 
temperature  rise  between  entrance  and  exit  of  the  heater. 

2.2.  Single  one  cover  solar  air  heater 


the  thermal  capacitance  effects  and  the  conduction  losses.  The 
energy  balance  for  each  component  of  the  collector  yields  the 
following: 


LUVCl  gldbi 


+ 


.1  6ia^  .  ugJ  —  ivig^g  -gp  I"  ivgug  nw\ 

hrgs(Tg  —  Ts)  +  hcgf(Tg  —  T  f)  +  hrpg(Tg  —  T  p) 


Absorber  plate  :  aprpS  = 
+  frcp  f(J  p  —  T  f)  +  hrpg(T 


MpCp¥  +  Vp 

p  —  Tg)  +  Ur(Tp  — 


d2Tp 
dx 2 
To) 


Airflow:  M  fC  f 


dT f  ,  GfCf  d2T f 


+ 


dt  w  dx2 
=  hCgf(T  f  —  Tg)  +  hcpf(Tp  —  T  f) 


(8) 


The  solar  radiation  data  and  ambient  temperature  have  been 
represented  by  Fourier  series  as  follows: 


5  =  237.17  +  388.7  cos(toi  -  3.272)  +  195.8  cos(2toi 
-  0.2757)  +  27.62  cos(3 tcc>\  -  3.8)  +  38.44  cos(4to! 

-3.44)  + 17.76  cos(5toi  -4.536)  (9) 

and 


Ta  =  14.354  +  6.1 1  cos(toi  -  3.858)  +  1 .88  cos(2toi 
-  0.509)  +  0.225  cos(3toi  -  1 .25)  +  0.288  cos(4toi 
-2.62)  + 0.146  cos(5to!  -3.92)  (10) 


with  oj\  =  (tt/12)  rd/h. 

The  initial  conditions  considered  by  these  authors  were: 
Tj(t  =  0 )  =  Tp(t  =  0)  =  Tg( t  =  0)  =  Ta(t  =  0 )  and  the  boundary  condi¬ 
tions  were  given  as  follow: 


97g 

dx 


x-0 


0, 


dTp 

dx 


=  0, 

x=L 


97g 

dx 


x=L 


0, 


dTp 

dx 


x=0 


0, 


(ID 


Eqs.  (6)-(8)  have  been  solved  numerically  using  explicit  finite 
difference  analysis.  The  results  showed  the  effect  of  various  design 
parameters  of  the  heater  on  its  performance. 

This  design  of  solar  air  heater  has  also  been  investigated  by  Ong 
[7].  Using  the  thermal  network,  the  following  equations  have  been 
written  for  a  single  one  cover  solar  air  heater  [7]: 


Cover  :  agS  +  hrpg(T P-Tg)  +  hcgf(T f  -  Tg) 

—  (hw  +  hrs)(Tg  —  Ta) 


(12) 


This  design  consists  of  single  channel  with  single  air  flow 
between  the  cover  and  absorber  plate  ( see  Fig.  2).  Many  authors  have 
investigated  the  thermal  performance  of  this  design  [7,9,13-18]. 

Garg  et  al.  [13]  developed  the  transient  heat  transfer  model  of  a 
single  one  cover  solar  air  heater.  This  transient  model  has  included 


The  wind  heat  transfer  hw  has  been  computed  using  (Al): 
Fluid  between  cover  and  absorber  :  hcpf(Tp  -T f) 

=  hcgf(T f-Tg)  +  2mcp  ^ ^  (13a) 
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with 

j  (■ T  f0  +  T  fi) 

Tf™= - 2 - 


(13b) 


Absorber  plate  :  rgapS  =  hcpf(Tp  -Tf)  +  hrpg(Tp  -  Tg) 

+  Ub(Tp-Ta)  (14) 

The  methodology  used  for  this  design  of  solar  air  heater  is 
similar  to  the  one  used  by  the  authors  to  study  a  bare  plate  solar 
heater. 

The  mathematical  model  proposed  by  Naphon  and  Kongtragool 
[17]  is  similar  to  the  one  of  Ong  [7]  instead  of  Eq.  (14)  which  is 
replaced  by 

dT 

mCp-^  =  hfg(Tf-Tg)  +  hfp(Tp-Tf)  (15) 

The  above  equations  proposed  by  Naphon  and  Kongtragool  [17] 
have  been  solved  by  the  explicit  method  of  the  finite  difference 
scheme.  The  results  obtained  from  the  model  have  been  compared 
with  the  experimental  data  of  Yeh  et  al.  [25].  It  appeared  that 
numerical  results  and  experimental  data  are  in  good  agreement. 

Mohamad  [15]  investigated  also  this  design  of  solar  air  heater. 
The  equations  proposed  by  this  author  are  similar  to  those  of 
Naphon  and  Kongtragool  [17],  instead  of  equation  for  the  cover 
which  did  not  contain  the  term  of  radiative  exchange  between  the 
glass  cover  and  an  ambient.  The  main  objective  of  this  author  was 
to  minimize  heat  losses  from  the  front  cover  of  the  collector  and  to 
maximize  heat  extraction  from  the  absorber. 

With  the  aim  to  determine  the  flow  channel  optimum  geometry 
of  this  single  cover  solar  air  heater  (Fig.  2),  Hegazy  [16]  proposed  a 
new  mathematical  model.  This  model  was  based  on  the  following 
assumptions:  (i)  the  flow  channel  was  hydraulically  smooth;  (ii) 
there  was  no  air  leakage  to  or  from  the  heater;  (iii)  heat  capacity 
effects  of  glass  cover,  enclosed  air,  absorber  and  bottom  plate  were 
negligible;  (iv)  the  temperatures  of  glass,  absorber  and  bottom 
plates  were  varied  only  along  the  x-direction  of  the  air  flow;  and  (v) 
the  solar  air  heater  consisted  of  a  number  elements  so  that  the 
temperatures  of  the  element  surfaces  are  uniform,  while  that  of  the 
air  stream  inside  varies  linearly  along  its  small  length.  With  these 
assumptions,  the  energy  balance  equations  have  been  written  as 

For  glass  cover  :  agS  +  hrpg(Tp  -  Tg)  +  hcgf(T fm  -  Tg) 

=  (hw  +  hrs)(Tg  —  Ta)  (16) 

where  Tfm  is  given  by  Eq.  (13b). 

For  air  stream  between  cover  and  absorber 

:  hcPf(TP  ~Tfm)=  hcgf(Tfm  -  Tg)  +  2GCP  (17) 

G  is  a  specific  mass  rate. 


At  absorber  plate  :  rgapS  =  hcpf(Tp  -  T fm)  +  hrpg(Tp  -  Tg) 

+  Ub(Tp-Ta)  (18) 

An  iterative  procedure  has  been  employed  by  the  author  to 
analyze  the  performance  of  this  solar  air  heater  in  terms  of  useful 
energy,  thermal  efficiency  and  air  temperature  rise. 

Hegazy  [26]  investigated  the  effect  of  variation  in  the  absorber 
width  on  both  the  thermal  and  hydraulic  performances  of  a  single 
cover  solar  air  heater  (Fig.  2).  In  addition  to  the  assumption  made 
by  the  author  [16],  it  was  assumed  that  heat  transfer  is  quasi¬ 
steady  and  one  dimensional,  and  that  turbulent  air  flow  inside 
small  element  temperatures  varies  linearly  along  its  small  length. 
Hence,  the  energy  balance  equations  of  this  design  of  solar  air 


collector  are  given  by  Eqs.  (16)  and  (18)  for  glass  cover  and 
absorber  plate.  For  the  air  stream,  the  energy  equation  has  written 
as 

hcpf(T p  -  T fm)Wdx  =  hcgf(T fi„  -  Tg)Wdx  +  mC„(T fo-Tfi)  (19) 

Ub  has  been  computed  using  the  same  equation  as  Hegazy  [16], 
in  which  W  =  W(x)  is  the  average  width  of  an  element  dx  at  a 
distance  x  from  the  inlet  and  G  =  m\Ap.  This  study  was  considered 
as  a  more  general  and  fundamental  case  when  variable  width 
absorbers  were  utilized,  instead  of  conventional  rectangular  plates, 
in  collecting  incident  solar  radiation.  The  results  indicated  the 
variable  width  collectors  with  shape  parameters  exhibit  perfor¬ 
mance  behavior. 

Aboul-Enein  et  al.  [27]  have  analyzed  this  solar  air  heater  device 
with  and  without  thermal  storage.  In  order  to  write  the  energy 
balance  equations,  these  authors  have  made  the  following 
assumptions:  (i)  the  heat  capacities  of  the  glass  cover,  absorbing 
plate  and  insulation  were  negligible;  (ii)  there  was  no  temperature 
gradient  across  the  thickness  of  the  glass  cover  and  the  storage 
material  has  an  average  temperature  at  time  t.  (iii)  the  system  was 
perfectly  insulated  and  there  was  no  air  leakage  from  the  collector; 
(iv)  no  stratification  existed  perpendicular  to  the  air  flow.  Without 
the  storage  material,  the  energy  balance  equations  of  the  glass  cover 
and  for  the  plate  are  similar  to  those  proposed  in  the  mathematical 
model  of  Hegazy  [16].  For  the  air  flow,  the  model  was  written  as 

hcpf(Tp  -  Tf)W  =  mCf^J-  +  dfp ,WC,£  +  hcfg (Tf  -  Tg)W 

+  Ut(T  f  —  Ta)W  (20) 

with  the  storage  material,  the  energy  balance  equation  for  the  glass 
cover  has  the  same  form  as  without  the  storage  material.  However, 
the  energy  balance  equations  for  the  absorber  plate  and  storage 
material  were  written  as  follow: 

For  the  absorber  plate  :  rgapS  =  hrpg(Tps  -  Tgs)  +  hcpf(Tps  -  T fS) 

+  “■  (7ps  —  Tst)  (21) 

ZP 

For  the  storage  material :  —  (Tps  -  Tst)Ap 

zp 

Jt 

=  MstCst-d±  +  Ub(Ab+As)(Tst-Ta)  (22) 

The  forced  convective  heat  hcga  is  given  in  Appendix  A. 

A  computer  program  was  prepared  for  the  solution  of  this 
mathematical  model.  An  optimization  process  for  a  flat  plate  solar 
air  heater  with  and  without  thermal  storage  was  carried  out. 
Comparisons  between  experimental  and  theoretical  results 
indicated  that  the  proposed  mathematical  model  could  be  used 
for  estimating  the  thermal  performance  of  flat  plate  solar  air 
heaters  with  reasonable  accuracy. 

Zhai  et  al.  [  1 8]  have  recently  investigated  a  single  cover  solar  air 
heater  integrated  with  the  building.  Based  on  some  assumptions, 
the  equations  governing  the  performance  of  the  solar  air  heater 
were  given  by  equations  similar  to  (12)-(14)  without  the  term  of 
the  radiative  exchange  between  cover  and  the  sky,  and  with 
Qi  =  mCp(Tfo  -  Tr)  and  Ub  =  [l/hCbattic  +  Sblkb]~\ 

Njomo  [14]  and  Njomo  and  Daguenet  [9]  analyzed  the  heat 
transfer  in  the  single  cover  solar  air  heater.  The  following  energy 
balance  equations  have  been  written  for  the  components  of  this 
design: 

For  combined  plastic-glass  transparent  cover: 


dT 

MgCPg^ 


g 


—  agS  +  hrpg(T  p  —  Tg)  +  hcgf(T  f  —  Tg) 


—  ( hcga  +  hrca)(Tg  —  Ta) 


(23) 
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_  .  „  _  dTp  mCpf  dTp 

For  air  flow:  PfefCpf-^  +  -Jf 

=  hcpf(Tp  -  Tf)  -  hcgf(Tf  -  Tg)  (24) 

dT 

For  absorber  plate  :  MPCPP  ^  —  (otx)S  hcp  f(T  p  Tf) 

-  hrpg(Tp  -  Tg)  -  (Ub  +  hrpa)(Tp  -  To)  (25) 

The  optical  efficiency  of  the  collector  is  given  by  [29]: 


Njomo  and  Daguenet  [9]  and  Njomo  [31  ]  have  also  investigated 
heat  transfer  in  this  design  (with  plexigas-glass  cover,  glass  cover- 
glass  cover  and  plastique-vitre).  These  authors  in  their  investiga¬ 
tion  have  considered  the  two  covers  as  a  combined  system  which  is 
characterized  by  one  heat  transfer  equation.  In  this  case  the  energy 
balance  equations  of  the  different  component  of  the  solar  air  heater 
design  are  given  by  Eqs.  (23)-(25)  without  the  radiative  exchange 
between  absorber  plate  and  sky.  (tot)  is  given  by  Eq.  (26)  with: 


_  Zg\  Xg2 

1  Pg\  Pg2 


(31a) 


(ra) 


X gOi  p 

1  -  (1  —  ot p) Pg 


(26) 


pg  pg\  +  i  _ 


Pg2  rgl 


Pgl  Pg2 


(31b) 


Tsky  was  given  by  the  correlation  proposed  by  Clark  [12].  This 
mathematical  model  has  been  used  by  Njomo  [30]  to  study  the  cost 
of  the  energy  extracted  from  these  collectors  as  a  function  of  their 
lifetime.  Numerical  calculations  have  been  performed  for  glass 
cover  and  plastic  cover  air  collectors  of  the  same  dimensions  under 
identical  conditions  of  insolation  and  mass  flow  rate.  The  results 
showed  that  a  significant  gap  exists  between  the  thermal 
performances  of  the  glass  cover  and  plastic  cover  air  collectors 
and  the  useful  energy  delivered  by  the  plastic  cover  air  collector  is 
cheaper  than  the  energy  extracted  from  the  glass  cover  air 
collector. 

2.3.  Single  two  covers  solar  air  heater 

This  collector  is  made  of  two  covers  and  a  single  channel  air 
flow  between  the  second  cover  and  absorber  plate,  and  insulation 
(Fig.  3).  This  model  has  been  studied  by  Njomo  and  Daguenet  [9], 
Mohamad  [15],  Naphon  and  Kongtragool  [17]  and  Njomo  [31]. 

Mohamad  [15]  investigated  heat  transfer  in  this  design.  Under 
steady  state  operating  conditions,  the  energy  balances  equations 
for  the  mentioned  collector  written  by  this  author  are 

For  the  top  cover  :  ag\S  =  hw(Tg i  -  Ta)  +  hCg-[g2  (Tgi  Tg2) 

hrg\g2 (Tgi  —  Tg2)  (27) 

For  second  glass  cover  :  ag2xg\S  =  hcg2 /(Tg2  -  T f) 

+  ^cgig2  (T g2  —  Tgi)  +  fi/-gig2  (Tg2  —  Tgi)  +  hrg2P(Tg2  —  Tp)  (28) 


oig  =  \  -  pg  -  xg  (31c) 

In  the  mathematical  model  studied  by  Njomo  and  Daguenet  [9] 
and  Njomo  [31],  the  heat  exchange  between  the  two  covers  was 
neglected. 

2.4.  Back-pass  solar  air  heater 

In  this  design,  the  absorber  plate  is  placed  directly  behind  the 
transparent  cover  with  a  layer  of  static  air  separating  it  from  the 
cover  (Fig.  4).  The  air  to  be  heated  flows  between  the  inner  surface 
of  the  absorber  and  the  layer  of  insulation  [3].  Heat  transfer  in  this 
design  has  been  investigated  by  Choudhury  and  Garg  [4],  Ong  [7], 
Hegazy  [26],  Garg  et  al.  [32],  Choudhury  et  al.  [33],  Al-Kamil  and 
Al-Ghareeb  [34],  Jannot  and  Coulibaly  [35],  Hegazy  [36]. 

Garg  et  al.  [32]  developed  a  simple  mathematical  model  to 
investigate  the  performance  of  the  system  under  different  design 
and  flow  conditions.  The  energy  balance  equations  for  different 
components  of  the  air  heater  were  written  as  follows: 

For  cover  :  hgs(Ts  -  Tg)  +  hrpg(J p-Tg)  =  Ut(Tg  -  Ta)  (32) 

For  absorber  plate  :  (ax )S  =  hcpf(T p  -  T fm)  +  hrps(T p  -  Ta) 

+  hrpb(T  p  —  Tb)  +  hrpg(T  p  —  Tg)  (33) 

For  the  flow  air  :  hcpf(Tp  -  T fm) 

=  hCbf(T fm  -  Tb)  +  mCa(T f0  -  Ta)  (34) 


For  air  stream  : 


mC, 


dT 


f 


dx 


^cg2/(Tg2  T  f)  +  hcpf(T  p  T  f) 

(29) 


Bottom  plate  :  hrpb(T P-Tb)  =  hcb f(Tb  -  T fm)  +  Ub(Tb  -  Ta) 

(35) 


For  absorber  plate  :  ctpXg\Xg2S 

—  hep /(Tp  —  Tf)  +  hrpg2(Tp  —  Tg2)  +  Ub(T  p  —  Ta)  (30) 

Naphon  and  Kongtragool  [17]  have  studied  this  design.  The 
mathematical  model  obtained  by  these  authors  is  similar  to  the 
model  of  Mohamad  [15]  with  change  in  the  equation  describing 
the  heat  transfer  in  the  upper  cover  in  which  the  radiative 
exchange  between  the  upper  cover  and  ambient  is  taken  into 
account. 


With  a  deal  to  compare  theoretical  parametric  analysis  of  solar 
air  heating  collector  with  and  without  packing  in  flow  passage 
above  the  rear  plate,  Choudhury  et  al.  [4]  investigated  this  design 
of  solar  air  heater  (Fig.  4).  The  steady  state  energy  balance 
equations  for  different  components  are  presented  below: 

For  cover  :  agS  +  hrpg(Tp  -  Tg)  +  hCpg(Tp  -  Tg) 

=  (hw  +  hrs)(Tg-Ta)  (36) 


hw 
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Fig.  3.  A  schematic  view  of  a  single  cover  solar  air  heater  [7,9,16,31]. 


Fig.  4.  A  schematic  view  of  a  single  cover  solar  air  heater  [7,16,26,31,32,36]. 
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For  absorber  plate  :  aprgS  =  hcpf(Tp  -Tf)  +  hcpg(Tp  -  Tg) 
+  hrpb(T  p  —  Tb)  +  hrpg(T  p  —  Tg) 


(37) 


For  the  flow  air  : 


T /)  =  hcbf(T f  -  Tb)  + 


mC  f  dT f 
~W~~dx 


(38) 


Bottom  plate  :  hcb f(T /  -  Tb)  +  hrpb(T P-Tb)  =  Uh(Th  -  Ta) 

(39) 

The  boundary  conditions  for  the  above  sets  of  equations  were 
Tj(x  =  0)  =  Tfit  Tj(x  =  L)  =  Tfo. 

In  the  mathematical  model  written  by  Ong  [7],  there  is  no  heat 
extracted  from  the  system  between  the  top  cover  and  the  stagnant 
air  cavity.  The  air  has  been  assumed  to  recirculate  within  the  space 
under  free  convection  in  an  inclined  wall  cavity.  Hence  the  energy 
balance  equations  obtained  for  this  design  are  the  same  as 
predicted  by  Choudhury  and  Garg  [4]  instead  of  Eq.  (38)  which  has 
been  replaced  by 

hcPf(TP  -  Tf)  =  hcbf(Tf  -  Tb)  +  {Tf  _  Tfi)  (40) 

In  the  investigation  of  Choudhury  et  al.  [4]  Eqs.  (36),  (37)  and 
(39)  should  be  considered  without  the  terms  of  natural  convection 
between  cover-absorber  plate  and  radiative  exchange  between 
cover-ambient.  Eq.  (37)  is  replaced  by 

1JWLdW  =  hpf^p-Td  +  h»f^-Tf)  (41) 

This  solar  air  heater  model  has  also  been  investigated 
experimentally  and  theoretically  by  Al-Kamil  and  Al-Ghareeb 
[34].  They  studied  experimentally  and  theoretically  the  effect  of 
solar  radiation  on  the  solar  collector.  The  balance  energy  equations 
based  on  the  subdivision  concept  were  developed  and  for  all  the 


components  of  the  solar  air  heater,  were  written  as 

Flowing  in  the  channel  :  hcpf(Tpi  -  T fi)Ai 

=  hcbf(T fi  -  Tbi)Aj  +  mCf(Tf(i_ i)  -  T p)  (42) 

Absorber  :  (ar)S  +  hcpf(Tpi  -Tf)  +  hnc(Tpi  -  Tgi) 

+  hrpb(Tpi  -  Tbi)  +  hrpg(Tpi  -  Tg,)  =  ppVpCp  T*  ~  (43) 

Glass  cover  :  agSFshFdAi  +  hrpg(Tpi  -  Tgi)Aj  +  hnc(Tpi  -  Tgj)Aj 
=  (hw  +  hrs)(Tgi  -  Ta)Ai  (44) 

For  bottom  plate  :  hcbf(Tfi-  TW)A  +  hba{Ta  -  Tw)/\, 

+  hrpb(T pi  ~  Tbi)Aj  =  0  (45) 


where  Fsh  and  Fd  are  shading  and  dust  factors.  The  balance  energy 
equations  were  solved  by  Gauss-Seidel  iterations  using  a 
computer  program  written  in  BASIC.  The  heat  transfer  coefficients 
were  determined  by  familiar  methods  available  in  the  literature. 
The  results  have  shown  a  closing  agreement  between  predicted 
and  measured  data  warrant  the  use  of  the  subdivision  method  to 
predict  some  useful  information  about  the  radiation  effect. 

Another  mathematical  model  for  this  design  of  solar  air  heater 
with  a  plastic  film  has  been  proposed  byjannot  and  Coulibaly  [35]. 
They  studied  theoretically  the  evolution  of  the  radiative  fluxes 
inside  the  collector  by  a  ray  tracing  method  applied  to  each  flux 
considering  the  different  and  successive  reflections  and  absorp¬ 
tions  occurring  on  absorber  plate  and  on  the  transparent  cover.  The 
thermal  analysis  will  result  in  the  radiative  balances  of  these  two 


components.  Assuming  that  the  collector  area  is  equal  to  1  m2,  the 
net  radiative  flux  received  by  these  two  components  are  stated  as 
For  absorber: 


<P 


Oi  ps'Egs 


s  + 


OlgiOi  pi 
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For  transparent  cover: 
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(46) 


(47) 


in  the  case  of  a  transparent  cover  opaque  to  IR,  rgi  =  0  and 
Pgi  =  1  -  agi.  Eqs.  (46)  and  (47)  become: 


0  _  OipsTgs 

1  PgsO  Oi  ps) 


S  +  <7 


P  g 


1  /  Oi  pi  1  /  Oigi  1 
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Tp  -  T4b 
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(48) 


and 


ip  g  —  Soigs  + 


Oigs  Egs  (1  —  Oi  ps) 
1  Pcs(l  Oi  ps) 


S  +  a 
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(49) 


The  results  showed  that  the  modifications  induced  by  taking 
into  account  the  partial  transparency  of  the  cover  towards  IR 
radiation  may  be  important. 

Another  investigation  has  been  published  by  Hegazy  [36].  In 
this  model,  Eq.  (16)  is  reconsidered  for  the  cover  balance  equation, 
with  Tfin{={Tfi  +  Tf0)l 2))  replaced  by  Tp  and  the  hg  becomes  hcpg. 
Energy  balance  on  the  bottom  plate  was  written  as  Eq.  (35).  Ub  is 
recalled  in  Table  1 .  For  the  others  components  of  the  design,  energy 
balance  equations  have  been  written  as 


For  absorber  plate  :  aprgS  =  hcpf(Tp  -  T fm) 

+  (hcpg  +  hrpg)(T  p  —  Tg)  +  hrpb(T  p  —  Tb) 


(50) 


For  the  air  flow  :  hcpf(Tp  -  T fm) 

=  hcbf(T fm  -  Tb)  +  2CCP (51) 

In  order  to  demonstrate  the  engineering  accuracy  of  the 
proposed  criterion,  the  author  has  determined  the  absorber  mean 
temperature  at  the  exit.  To  accommodate  this  necessity,  energy 
balances  (35),  (50)  and  (51)  have  been  applied  progressively  to 
each  element  and  solved  iteratively.  It  was  concluded  that  the 
channel  ratio  (D/L)  plays  a  decisive  role  in  determining  the  rate  of 
useful  heat  gain. 

To  determine  the  flow  channel  optimum  geometry  and  study 
the  effect  of  the  variation  in  the  absorber  width  on  both  the 
thermal  and  hydraulic  performance  of  this  design,  two  new 
mathematical  models  have  been  published  by  Hegazy  [17,26].  In 
these  new  models,  all  energy  balance  equations  proposed  by  this 
author  in  1999  are  reconsidered  instead  of  the  equation  of  the  air 


Table  1 

A  review  of  heat  transfer  coefficient  used  in  different  models. 


References  Solar  air  heaters  systems 
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Aboul-Enein  et  al.  [27] 


Njomo  [14]  and  Njomo  and  Daguenet  [9] 
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Parallel-pass  solar  air  heater:  Type  2 


Heat  transfer  coefficients 


The  radiative  and  natural  convective  coefficients  were  summarized  in  Duffie  and  Beckmann  [5] 

hPf  and  hbf[ 6]:  h/=Nuoc(l  +MD/L)/c/D  with  Nu^  =  0.0182Re°8Pr°4  and  M  =  14.31og10(N)  -  7.9.  IV  =  L/D  if  0  <  L/D  <  60 

and  N  =  60  if  L/D  >  60. 

hcpa  =  0.86 Rel/2Pr1/3 (kp/Lp),  with  Lp  =  4WL/2W  +  2 L  for  rectangular  plates  and  with  2  x  104  <  Re  <  6  x  104  [10]. 

Duffie  and  Beckmann  [5]  suggested  that  the  above  equation  is  valid  for  Re  >  10s. 
hcpf  and  hcbf( see  [11]):  for  the  laminar  flow  (Re  <  2100):  Nu  =  hcDh/I< f  =  4.9+ 

(0.0606(RePrDh/L)12/l  +  0.0909(RePrDh/L)°  7Pr° 17)  for  turbulent  flow  (Re  >  2100)  by  [5]:  Nu  =  hcDh/I<f  =  0.0158 Re08 
with  Dh  =  4WD/2W+2D 

hr pa  =  (j£p(Tp  -  T^ky)/(Tp  -  Ta ),  where  Tsky  =  TJ0.787  +  0.0028(Td  -  273 )]025  [12]  and  Td  is  the  dew  point  temperature 
of  the  ambient  air;  it  is  related  to  its  water  vapor  pressure  emb  (in  millibars)  and  has  been  calculated  from: 
emb  =  1013.25/760  exp(20.519  -  5179.25/Td),  hrbp  =  cr(Tp  +  Tb)(T2p  +  T2)/( 2/ep)  -  1 

The  important  dimensions  and  the  heat  transfer  parameters  used  were  given  in  McAdams  [19],  Kreith  and  Kreider  [20] 
The  radiation  coefficient  between  parallel  plates  set  i-j  is  given  as:  hnJ-  =  er(T;  +  Tj)(Tf  +  P2)/(  1  /eft  +  (1  /eft  -  1 
The  forced  convection  heat  transfer  coefficient  between  parallel  plates  is  given:  for  laminar  flow  (Re  <  2300) 

[21]  hDh//c  =  5.4  +  0.00190(RePrDh/L)171/l  +  0.00563(RePrDh/L)171,  for  transition  flow 

region  (2300  <  Re  <  6000)  [22]:  ( hdh/k )  =  0.116(Re2/3  -  125)Pr1/3[l  +  (Dh/L)2/3](/i/pw)° 44  and  for  turbulent  flow 
(Re  >  6000)  [23]  (hDh/k)  =  0.01 8  Re0  08  Pr04,  the  heat  transfer  coefficient  at  the  entrance  region  is  computed  with 
the  same  relation  as  hpf( see  [4]),  hrs  =  <JSg(Tg  -  T4)/(Pg  -  Ta)  ^ 

The  convective  heat  transfer  coefficient  is  computed  using  [24]:  Nu  =  (hDh/fc)  =  0.333  Re08  Pr1  3  with 
Dh  =  4WB/2W  +  2B 

hcpf  and  hcgf  are  computed  with  [6]:  h  =  fc/Dh{0.0158Re°8  +  (0.00181ReDh  +  2.92)exp(-0.0379x/Dh)} 
where  ReDh  =  2M/(1  +  D/W )  Dh  is  hydraulic  diameter  and  M  is  a  mass  flow  parameter.  Ub  =  [1  +  2 (8in  +  D 
+  Sg/WWinlkin  +  1  lhw]~\  8in  and  8g  are  the  thickness  of  insulation  and  glass  cover. 

The  radiative  heat  transfer  coefficient  (see  Ong  [7]) 

hcpf  and  h Cfh  were  computed  using  [28]:  Nu  =  hcpgDhlk  =  1  +  1.44[1  -  1708/Ra  cos(/3)]+  +  [1  -  1708(Sin(1.8/3))1(5/Ra 
Cos(ft)]  +  [(Ra  cos(/3)/5830)1/3  -  1  ]+  where  [  ]+  means  that  only  positive  values  are  taken  into  account,  fi  is 
a  collector  tilt  angle. 

hrpg  (see  Ong  [7])  and  hrgs  =  crsg(T j  +  T2kys)(Tg  -  Tsky ) 
hcgp ,  hrpa,  hcpf  (see  Choudhury  et  al.  [4]) 

hrij  =  ci8p((  1  -  Pg)Tp  -  /lgT4)/(l  -  ppPg)(Tp  -  Tg)  and  hrpa  =  aspTir(T4p  -  T%y)/(TP  -  Ta) 
hrgia  =  cr(Tgi  +Ta)(T21  +Ta)/((l/figi)  -  !).  hrpg( see  [7])  with  i  =  p  and  j  =  g 
The  radiative  heat  transfer  coefficients  (see  Ong  [7]) 

The  convective  heat  transfer  coefficient  from  the  upper  plate  to  the  flowing  air  and  the  pressure  drop  [37]: 

Nu  =  2 bh/k  =  Nu0  +  P(b/L),  A P  =  [/0  +  7(ib/L)(m2/p)(L/fa)3],  for  Re  <  2550  (laminar  flow)  Nu0  =  5.385; 

/3  =  0.0148Re;  /0  =  24/Re;  y=  0.9;  for  2550  <  Re  <  104  (transitional  flow)  Nu0  =  4.4  x  10  4Re:  2; 

ft  =  9.37  Re0  471;  /0  =  0.0094,  y  =  2.92  Re  015 ;  for  104  <  Re  <  105  (early  turbulent  flow)  Nu0  =  0.03Re°  74; 

ft  =  0.788 Re0471;  /0  =  0.059 Re'02;  y=  0.73. 

hcpf  and  hcbf( see  [17])  other  heat  transfer  coefficients  (see  [5]) 

ub=  [1  +  (2/WXSi„  +  D  +  S  -  SgWtolkto  +  1  //!„]-’,  hcpg  (see  [27]) 

hcpf  and  hcbf  is  correlated  by  (see  Biehl  [38]):  h  =  /</Dh{0.015Re^  +  (0.001 81ReDh  +  2.92)exp(-0.03795x/Dh)} 

The  forced  convection  heat  transfer  coefficient  between  air  and  plate  or  between  air  and  cover  is  obtained 
by  using  the  equation  [24]:  h  =  0.0293  Re08kf/D.  The  natural  convection  heat  transfer  coefficient  (see  Wong  [41]) 

The  radiative  heat  transfer  coefficients  and  the  forced  convection  heat  transfer  between  parallel  plates  for  laminar 
flow,  transition  flow  region  and  turbulent  flow  (see  Ong  [7]).  For  9500  <  Re  <  22000,  situations  involving  a 
large  property  variation  [42]  Nu  =  ( hDh/k )  =  0.027  Re08Pft/3(ii/pw)0A4 

hcpf i,  hcgfi,  hcbf2,  hcpf2  were  calculated  using  the  following  equation:  h  =  (k/Ddj){0.015Re°8  +  (0.001 81Rej 
+2.92)exp(-0.03795x/Ddj)},  Ub  =  [1  +  2 (8in  +  D^+D2  +  8gIW)][8inlkin  +  l/hw]_1  W  was  replaced  by  the 
average  width  VV(x)  in  the  second  model  [26] 

Radiative  heat  transfer  coefficients  (see  [7-9]),  free  convective  heat  transfer  coefficient 
(see  [16]),  Uej  =  (UA)edgeIAc  =  {(klx)Sj}lW 

{klgia  =  Ujl  -  (hrpgl  +  hcpf1)  |  or  Ug^a  =  (hw  +  hrg 2fl)  -  (hCgig2  +  ^rg\g2) 

where  hCg\g2  —  1.25(7 fm\  1+2)  >  hyPg\  ~  AcfT)ml/W  +  8gl  ~  hrpb  «  4(jT}m2/K'  +  8b'  ~  ^ 

hrg^g2,  hw  (see  Ong  [7]  and  (Al)),  hrg2a  =  cr£g(Tg2m  +  Ta)(T22m  +  T2) 

Ug 1  =  1  +  hck f2  +  Ifikg1  and  Ub  ~  kslc  1 
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Matrawy  [49] 
Ammari  [51] 


Solar  air  heater  with  slats 


Wijeysundera  et  al.  [57] 

Verma  et  al.  [59] 

Two-pass  solar  air  heaters 

Naphon  and  Kongtragool  [17] 

Zhai  et  al.  [18] 

Yamali  and  Solmus  [63] 


Ho  et  al.  [66] 

Double-pass  flat-plate  solar 
air  heater  with  recycle 

Naphon  [71] 

Double-pass  flat  plate  solar  air 
heater  with  longitudinal  fins 

Chandra  et  al.  [72] 

Triple  pass  solar  air  heaters 

Ho  et  al.  [76] 

Multi-pass  flat-plate  solar  air 
heaters  with  external  recycle 

Tchinda  [78] 

Solar  air  heater  with  a 
compound  parabolic  concentrator 

Karim  and  Hawlader  [83] 

v-Groove  solar  air  collector 

Bashria  et  al.  [84-86] 


Choudhury  and  Garg  [88] 

Single  pass  solar  air  heater 
with  packed,  type  1 

Choudhury  and  Garg  [88] 

El-Sebaii  et  al.  [99] 

Ramadan  et  al.  [96] 

Single  pass  solar  air  heater 
with  packed,  type  2 

Double-pass  solar  air  heater 
with  packed 

The  forced  convection  heat-transfer  coefficient  between  two  flats  plates:  for  laminar  flow  in  a 
short  conduct  [21]:  Nu/  =  Dhhj/k  =  4.4  +  0.00398(0.7RejPrDh/L)166/l  +  0.0114(Re,PrDh/L)112  (j  =  1,  2), 
while  for  turbulent  flow  [19,46]  data:NUj  =  ( Dhhj/k )  =  0.018Re°os[l  +  (De/L)0  7] 

Ut=  [l/(hnc+  hrlc)+  l/(hw+  hrca)]-1 


hrca  =  crsc(Tc  +  TS)(TS2  +  T2)(Tc  -  Ts/Tc  -  Ta),  hrpg  and  hrga  (see  [7]),  Natural  convection  (see  [27]),  inner 
surface  convective  coefficients  (see  [7-9]) 

U21  =  [(l/ftest /)  +  (l/hcg2f)]"1,  hrij  =  ohij (rf  +  T] )  (T,  +  Tj),  hdj  =  Cu( T,  +  Tj). 

The  convective  heat  transfer  coefficient  between  enclosed  spaced  (see  [27]) 

The  forced  convective  heat  transfer  coefficient  (see  [8,9]) 

The  heat  transfer  coefficient  (Niles  et  al.  [24]):  Nu  =  hD/k  =  0.0333  Re08Pr1/3. 

hg\ a  =  convection  heat  transfer  from  glazing  due  to  wind  (see  Eq.  (Al))  +  radiation  heat  transfer 

coefficient  from  glass  cover  to  sky  referred  to  the  ambient  air  temperature(see  [8,9]). 

The  global  heat  transfer  coefficient  from  insulation  plate  to  the  attic:  Ub  =  [l/hCbattic +  Ab//<fa] 
frisky  (see  [8,9]),  hrg\g2  (see  [7]),  hcg\g2  (see  [28]),  hcg\f\  (see  [8,9])  for  laminar  flow  and  (see  [64]): 
Niw,  =(/;,/ 8)  {Rep  -  1000)Pr/l  +  12. 7(/ /i/8)°  5  (Pr0,67  -  1)  with  fn  =  ^0.79  In  Re,, -  1.64. 
h1  =  h\,  h2  =  h'2.  The  forced  convection  heat  transfer  between  two  plates  (see  Kays  [46]),  for  a 
short  conduit  h}  and  h2  (see  [19]):Nuj  =  hjDe/k  =  0.0158 R° 8  1  +  (De/L)07  (j  =  1,2),  for  the 
laminar  flow  (see  Heaton  et  al.  [21]),  hT( see  [67,68]  and  Eq.  (A4)). 

1/hfi  =  1/hfiis  +  1/^2  +  h/ks  VUgis  =  l/hrg2s  +  hw  +  l/hCgig2  +  ^rgig2  or  l/L/g/s  =  1 IUT+  +  h\ 

where  hcgUg2  =  1.25 (Tglm  -  Tg2m)a25  (see  [67]),  hrpgl  =  4crT3m/[(l/ep)  +  (l/ep)  -  1] 

hrpR  =  4crT3m/[(l /ep)  +  (1  /er)  -  1],  those  between  two  glass  covers  (see  Ong  [7])  and  from 

cover  2  to  ambient:  hrg2s  =  (J8g(T^2  +  T2(y)(Tg2  +  Tsky ). 

The  convective  heat  transfer  coefficient  between  channel  (see  [7,17]). 


The  radiative  heat  transfer  coefficient  between  two  surfaces  i  and  j  (see  Ong  [7])  and  the 
convective  heat  transfer  coefficient  (See  [73]):  Nu  =  h / D/k  =  0.036  pr08pr08 
h1  =  h\ ,  h2  =  h2,  h3  =  h'3,  andh4  =  h'4 , .  Other  coefficients  (see  Ho  et  al.  [66]) 

hn  =  cr(T}  +  T2)(Tc  +  r„)/l/e„  +  Av/Ac(l/ec  -  1)AC/AP,  hr2  =  G8C{T2C  +  T2)(TC  +  TS)AC/AP,  hv/c  =  (3.25 
+0.0085T„  -  Tc/2Dv)Av/Ap  (see  [79])  with  Dv  =  2lv(ef+  epv)llv  +  ef+  epv  Uf( see  [5,11,81]) 

The  convective  heat  transfer  coefficient  within  the  airflow  duct  (see  [61]):  Nu  =  Nu0  +  n/3b/L,  For  laminar 

flow  (Re  <  2800):  Nu0  =  2.821,  /3  =  0.126  Re;  For  transition  flow  (2800  <  Re  <  105):  Nu0  =  1.9  x  10_6Re1  79,  fi  =  225; 

For  turbulent  flow  (104  <  Re  <  105):  Nu0  =  0.030210  6Re0  74,  /3  =  0.242  Re0  74,  the  top  loss  coefficient  Ut, 

(see  Klein  [67]) 

The  radiative  coefficient  from  absorber  to  the  cover  [52,87]: 
hr  =  <j(J2p  +  Tj)(TP  +  Tg)/\/Spf  +  l/eg  -  1  with  ep/  =  2ep/l  -  ep; 

The  convective  heat  transfer  (see  [7]),  with  Dh  =  2Hysin(0/2)/l  +  sin(0/2) 

hcpf,  hcbf  and  hcpl/,  (see  [89]):hcp/  =  hcbf  =  0.2Re°pfPr1/3k fD~\  hcpl  f  =  [l/hpl  f  +  DP/Wkpl]~\ 

with  hpl  f  =  0.255Re°f  e-1.  For  the  other  heat  transfer  coefficients  (see  [5]). 

For  all  the  heat  transfer  coefficients  (see  [88,90,91,95]) 

For  the  various  heat  transfer  coefficients  (see  [97]) 

hcg\  f-iNup-ikfi/De  =  0.2Re°pfPrl/3kfi/De  (see  [98]),  with  Repl  =  mfDe/fiAp i,  De  =  (2/3)[eDpi/(l  -  e)],  e  is  porosity  of  the 

packed  bed  (e  <  l):e  =  ( Vchan  -  Vp\)/Vchan,  Vchan  and  Vpl  are  the  total  volume  of  the  channel  and  the  volume  of  the 

packed  bed  material.  Dpi  =  6VR  3 /rm  (see  [97])  Vs  =  total  volume  of  n  particles  selected 

randomly.  Apl  =  6(1  -  e)/Dpl  (see  [97]).  hp/1  =  begin  (see  [43]),  hcpl/1  =  [l/hpl/  +  De/DplSp-l]~'1 

(see  [90])  with  Npl  f  =  hpl  fDe/kfi  =  (0.255 /s)Rep/3Pr^/3,  5pl  is  a  constant  and  depends  on  the 

packed  bed  shape. 
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Fig.  5.  (a)  A  schematic  view  of  a  single  cover  solar  air  heater  [7,16,26,32,36].  (b)  A 
schematic  view  of  a  single  cover  solar  air  heater  [45]. 


flow  which  has  been  replaced  respectively  by  Eq.  (17)  and: 
hP(Tp  -  Tfm)  =  hcbf(Tfm  -  Tb)  +  2GCP  (^’/~^)  (52) 

with  Re  =  4m/juP(x )  and  Dh  =  4 W(x)D/P(x). 

2.5.  Parallel-pass  solar  air  heater 
2.5.1.  Type  l 

The  solar  air  heater  configuration  under  consideration  is 
shown  in  Fig.  5a.  It  consists  of  one  cover,  double  channel 
design  with  double  air  flows  between  cover  and  absorber 
plate  and  between  absorber  and  bottom  plates  and  with 
insulation  provided.  Many  papers  have  investigated  this  design 
[7,16,17,26,39,40,43]. 

Jha  et  al.  [43]  obtained  the  energy  balance  equations  in  partial 
differential  forms,  which  govern  the  behavior  of  the  system,  by 
considering  energy  conversation  at  each  component  of  the  system 
separately.  For  the  solar  heater  shown  in  Fig.  5a,  the  following 
equations  were  obtained: 

3T 

At  the  cover  :  MgCg-^f  =  agS  +  hrpg(T p  -  Tg )  +  hc/lg(T fl  -  Tg) 
—  hCgw(Tg  —  Tw)  —  hrga(Tg  —  T  a)  (53) 


The  boundary  and  initial  conditions  are  given  by  the  following 
equations: 


dTp 

dx 


x=0 


0, 


dTp 

dx 


x=L 


o, 


dTb 

dx 


x=L 


o 


(58a) 


Tfl  (x  =  0)  =  Tfl,  Tf2(x  =  0)  =  Tfi 


(58b) 


Ong  [7]  obtained  from  the  thermal  network  the  following  heat 
balance  equations: 

For  cover  :  agS  +  hrpg(T P-Tg)  +  hcgf ,  (T ;i  -  Tg) 

=  {K  +  hrs)(T g  -  T a)  (59) 

For  air  flow  between  cover  and  the  absorber  plate  : 

hcPf i  (TP  —  Tfj)  =  hcg/1  (T ;i  -  Tg)  +  2 mtCf  Tf'w^f'  (60) 

For  absorber  plate  :  aprgS  =  hcpf2(Tp  -T f2)  +  /irpg(Tp  -  Tg) 

+  hrpb(T  p  —  Tb)  +  hcpfi(T  p  —  T  fi)  (61) 

For  air  flow  between  absorber  plate  and  the  bottom  plate  : 

hCpf2(Tp  -Tf2)  =  hcbf2(Tf2  -  Tb)  +  2  m2CfTf2w*if2  (62) 

With  the  aim  to  determine  the  channel  optimum  geometry  of  the 
flat  plate  and  to  study  the  effect  of  the  variation  in  the  absorber 
width  on  both  the  thermal  and  hydraulic  performances  of  this  solar 
air  heater  (Fig.  5a),  Hegazy  [16,26]  has  written  two  mathematical 
models.  In  the  two  models,  the  energy  balance  equations  on  the 
glass  cover,  the  absorber  plate  and  the  bottom  plate  are  given  as 

otgS  +  hg(T  fm- 1  —  Tg)  +  hrpg(T  p  —  Tg)  =  (hw  +  hrgS)(Tg  —  Ta)  (63) 

otpXgS  =  hcpfi(Tp  —  T  /mi)  +  hrpg(T  p  —  Tg)  +  hcpf2(Tp  —  T  fm2) 

+  hrpg(T  p  —  Tb)  (64) 

d-cb/2 (T fm2  ~  Tb)  +  hrpb(Tp  -  Tb)  =  Ub(Tb  -  Ta)  (65) 

The  energy  balance  equation  of  the  air  flow  in  the  upper  and  the 
lower  channel  are 

For  the  first  model  :  h  P/1  ( TP  -  Tfinl ) 

=  hg(jm  -  Tg)  +  2G1CP  r^/~^'71  (66) 

hpf2(.Tp  -Tfm2)  =  hb f2 (T fm2  -  Tb)  +  2G2Cp/^7i£  (67) 


For  air  flowing  between  cover  and  plate  absorber  : 

Mf'' cf  ■ =  z^r£ {Tp  ~T^~ ~ Tg) 

(54) 


For  Plate  1  :  MbCp 


9T, 


-  Oi  r)  T  crS  k  n  d  ; 


92T, 


hrpg(T  p  —  Tg) 


- PLgu  -  P  P~dx2  “rpgK1  p 

-  hcpf2(.T p  —  T f2)  —  hrpb(T p  -  Tb)  —  hcpf](Tp  —  T fi) 

For  air  flowing  between  absorber  plate  and  bottom  plate 

9T/2  -C2Cf  dTf2 


(55) 


W  dx 


+  hCp/2(T p  —  T f2)  —  hcbf2(Tb  -  T f2) 


(56) 


For  the  second  model 

T frn\  -  Tift 


hpft(Tp  —  Tfini)  —  hg(T fin i  ~  Tg) 


+  2Gi  Cp 


AP/A(x)) 


^  p/2  (Tp  —  T  fin 2)  —  hbf2(Tfin 2  —  Tb)  +  2 G2Cp 


T  f m2  -  T if2 

Ja p/m ) 


(68) 

(69) 


with  Re  =  4 m//xP  and  Dh  =  4W/P  for  model  1  [16],  Re  =  4 m//xP/(x) 
and  Dh  =  4 W(x)/Pj(x)  for  the  second  model  [26]. 

To  study  the  performance  of  this  solar  design  (Fig.  5a),  Naphon 
and  Kongtragool  [17]  proposed  another  mathematical  model.  This 
model  is  characterized  by  the  following  equations: 

For  the  cover  :  ag^S  =  hw(Tg i  -  Ta)  +  hc/lgl (Tgl  -  T n) 

+  h-rgi  p(Tg\  —  T  p)  +  hrag  i  (Tgi  —  Ta)  (70) 


For  bottom  plate  :  MbCb 


dTb 

dt 


=  -kh8 


d2T 


bdb  ~^2~  +  hrpb(T p  -  Tb) 


hCbf2{Tb  -  T f2)  -  hb(Tb  -  Tr) 


(57) 


For  the  air  flow  between  upper  cover  and  absorber  :  mC f 

—  ^cgi  ft  ( Tgi  —  T fi )  +  hcpf  \  (Tp  —  T/i ) 


dTn 

dx 

(71) 
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For  the  absorber  plate  OL  p  7 g\  S 

=  /I  p(^p  —  T'/i )  +  /2 p (T p  ~~  T  f 2)  +  hrgip(T p  —  Tg\  ) 

+  hrpb(T  p  —  Tb) 

For  the  air  flow  in  lower  channel  :  mCy-^- 

=  he f2p(T p  -  T f2)  +  hCf2b(Tb  -  T f2) 


(72) 

(73) 


Sp  are  given  by  the  following: 


otgRKj 


Ho 

3600 Nds 


f /?  a"  +  fi' 


n\  H0 
Nj  3600 Nds 


(82) 

(83) 


For  the  bottom  plate  :  hcf2h(Jb  -  T f2)  +  hrpb(Jb  -Tp) 

+  Uh(Th-Ta)  =  0  (74) 

Forson  et  al.  [43]  have  investigated  theoretically  and  experi¬ 
mentally  heat  transfer  in  the  design  shown  in  Fig.  5a.  Making  a 
general  assumption,  Forson  et  al.  [43]  formulated  the  equations 
governing  the  performance  of  the  system  by  coupling  the  energy 
equations  of  the  components  of  the  solar  air  heater  with  those  for 
the  useful  heat  extracted  in  the  two  channels.  The  energy  balance 
equations  for  various  components  of  the  system  are  given  below: 


For  the  cover  :  Sg  =  hw(Tgl  -  Ta)  +  hrcs(Tg  -  Ts) 

+  hrpg(Tg  —  T  p)  +  hcgf-[(Tg  —  T  (75) 

For  the  air  flow  between  upper  cover  and  absorber: 

m'Cfibr= W/Icp/1  {Tp  ~Tf')+ whcgf' {Tg  ~ T/1 } 

+  Uegsl(Ta-Tf])  (76) 

For  the  absorber  plate: 

5p  =  hcp f2(T p  —  T f2)  +  hcpfi  (T p  —  T fi)  +  hrpb(Tp  —  Tb) 

+  hrpg(Jp-Tg)  (77) 

For  the  air  flow  in  lower  channel: 

dT 

m2Cf-3j±  =  Whcf2b(Tb  -Tf2)  +  Whcpf2(Tb  -Tf2) 

+  Ues2(Ta-Tf2)  (78) 


with 


24  x  3600 

71 


(l  +  0.033 Cos 


/360 n\\ 


+  +  Sin(0)Sin(S) 


Cos(@)Sin(S)sincos 

(84) 


A  FORTRAN  computer  program  was  written  and  used  to  solve 
the  above  set  of  equations.  In  the  solution  procedure  it  was 
assumed  that  the  air  heater  is  subject  to  a  constant  insulation  and 
hence,  the  plate  temperature  is  constant.  It  was  shown  that 
significant  improvement  in  solar  air  heater  performance  can  be 
obtained  with  an  appropriate  choice  of  the  collector  parameters 
and  the  top  to  bottom  channel  depth  ratio  of  the  two  ducts.  The 
mass  flow  rate  was  shown  to  be  the  dominant  factor  in 
determining  the  overall  efficiency  of  the  heater. 


2.5.2.  Type  2 

The  energy  flow  diagram  of  such  a  device  is  presented  in  Fig.  5b. 
At  seen  in  Fig.  5b,  two  air  streams  of  different  flow  rates  but  with 
total  flow  rate  fixed,  are  flowing  steadily  and  simultaneously 
through  two  separate  channels  (above  and  below  the  absorbing 
plate).  With  the  presence  of  the  second  cover,  this  device  is 
different  to  the  one  of  Fig.  5a.  Yel  et  al.  [45]  investigated  heat 
transfer  in  this  design.  In  the  mathematical  model  proposed  by 
these  authors,  the  following  assumptions  have  been  considered: 
the  temperatures  of  absorbing  plate,  bottom  plate  and  bulk  fluids 
were  functions  of  the  flow  direction  only,  and  both  the  covers  and 
fluid  do  not  absorb  radiant  energy.  The  energy  balances  within  the 
differential  length  dx,  were  given  as 

For  inner  cover  :  /irpci  (Tp  -  Tgl )  +  hcgl  n  (T n  -  Tg  1 ) 

=  Ugla(Tgl-Ta)  (85) 


For  the  bottom  plate: 

hCf2b(Tb  -  T f 2)  +  hrpb(Tb  -  Tp)  +  Ub(Tb  —  Ta)  =  0 


Absorbing  plate  :  aprp  =  Ut(Tp  -  Ta)  +  UB(TP  -  Ta) 

(79)  +  hCpFi(Tp  -  T +  hCpF2(Tp  -  T f2)  (86) 


The  pressure  drop  is  given  by  Achenbach  and  Cole  [44]: 

Dt  =  233.9 hpa  ~  Ta)  (80a) 

with 

Tfm=^[Ta-0.65(Ta-Tfo)]  (80b) 

and 


For  the  bottom  plate  :  hcb f2(T f2  -  Tb)  +  hrpb(T p  -  Tb) 

=  Ub(Tb  -  Ta)  (87) 

Jt 

For  fluid  1  :  rmC f-£-=  Whcpn(Tp -T n) 

+  Whcgfi  (Tg  —  T fi)  +  hc/igi  (T f\  —  Tg\)  (88) 

JT 

For  fluid  2  :  (1  -r)mCf—^- 

=  WhCPF2(Tp  -  T/2)  +  Whcbf2(Tb  -Tf2)  (89) 


h  =  Lsin(/3) 


(80c) 


2.6.  Solar  air  heater  with  slats 


The  mass  flow  rate  of  air  in  each  of  the  two  channels  is  given  by 


mj  =  wPjSjCvj 


2D 


hj 


\  M 


233.9hpa 


fo 


J  fo  + 


(81) 


Dh j  is  the  hydraulic  mean  depth  of  the  duct  and  CVJ  the 
coefficient  of  the  velocity.  The  edge  heat  losses  are  estimated  by 
assuming  one  dimensional  sideways  heat  flow  around  the 
perimeter  of  the  collector  system.  In  Eqs.  (75)  and  (77)  Sg  and 


Yeh  et  al.  [47]  have  presented  a  mathematical  model  for  double 
flow  solar  air  heaters  with  fins  attached.  Fig.  6a  shows  the  design  of 
a  solar  air  heater  with  fins  attached  in  which  the  absorbing  plate 
divides  the  air  conduit  into  two  parts,  channel  1  and  2.  The  energy 
flow  diagram  of  such  a  device  is  also  presented.  As  seen,  two  air 
flow  of  different  mass  flow  rate  but  with  total  rate  fixed,  flow 
steadily  and  simultaneously  through  two  such  separated  channels 
of  the  same  size  for  heating.  The  method  for  theoretical  prediction 
of  collector  efficiencies  as  well  as  the  experimental  procedure  was 
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Fig.  6.  (a)  Schematic  view  studied  by  Yeh  et  al.  [47].  (b)  Schematic  view  studied  by 
Ammari  [51].  (c)  Schematic  view  studied  by  Kabeel  and  Mecarik  [50]. 


similar  to  that  presented  by  the  same  authors  [45],  except  that  the 
allowance  was  for  the  fins  to  be  attached  to  the  upper  and  lower 
surfaces  of  the  absorbing  plate.  It  was  shown  that  providing  fins 
attached  on  the  collector,  will  improve  the  collector  efficiency  and 
that  the  application  of  the  concept  of  double  flow  in  the  design  of  a 
solar  air  heater  with  fins  attached  was  technically  and  economic¬ 
ally  feasible. 

Al-Nimr  and  Damseh  [48],  Matrawy  [49],  Kabeel  and  Mecarik 
[50]  and  Ammari  [51]  have  investigated  a  solar  air  heater  system 
that  has  metal  slats  connecting  the  absorber  plate  to  the  bottom 
plate  as  a  development  to  enhance  the  thermal  performance  of  the 
solar  heater  (Fig.  6b). 

With  the  aim  to  study  the  dynamic  performance  of  this  design 
of  solar  air  heater,  Al-Nimr  and  Damseh  [48]  developed  a 
mathematical  model  on  the  base  of  the  following  assumptions: 
(i)  there  are  two  separate  phases,  fluid  and  solid,  which  are  not  in 
thermal  equilibrium.  This  was  the  case  for  low  thermal  con¬ 
ductivity  fluid  and  low  volumetric  heat  transfer  coefficient 
between  the  fluid  and  solid  domains;  (ii)  the  fluid  was  one 
dimensional  plug  flow;  (iii)  the  fluid  was  laminar  and  Newtonian; 
(iv)  the  properties  of  the  fluid  and  solid  domains  were  constant  and 
uniform;  (v)  thermal  and  mass  diffusion  were  assumed  to  be 
negligible;  (vi)  transient  response  resulted  from  sudden  change  in 
the  intensity  of  the  incident  solar  radiation  and/or  inlet  fluid 
temperature.  This  model  is  described  by  the  following  equations: 

9ft  =  WSLpfAf  _WLpfAfUu  _  bLpjAfU ,  _ 

9r  T^psAsCstn  f  psAsCsm  f  psAsCsm  f  >’ 


dr 


+ 


Wf 

dx 


dLUs  bLUi 

Csm  f  f  +  Csm  f 


(Os -Of) 


(91) 


where 

p  ( Tf  ~  Too) 

Of  = - f - . 

i  rv) 


X  = 


ft  = 


(Ts  ~  Tee) 


Oo 


(T f(t,0)  -  Too) 


OO 


oo 


tm 


f 


PfAf/L 


(92) 


with  the  following  initial  and  boundary  conditions:  0/ 0,x)  = 
0s(O,x)  =  0  and  0j{ r,0)  =  60(t). 

The  authors  have  outlined  analytical  solutions  for  the  fluid  and 
solid  domains.  The  effects  of  different  design  parameters  on  the 
thermal  performance  of  the  heater  were  analyzed  and  the  validity 
of  the  theoretical  model  was  verified  experimentally.  It  was  found 
that  both  theoretical  and  experimental  results  were  in  good 
agreement. 

Matrawy  [49]  proposed  another  mathematical  model  of  this 
solar  air  heater  design.  The  model  was  based  on  the  following 
assumptions:  (i)  uniform  temperature  distribution  across  the 
width  of  absorber  plate  and  of  back  plate;  (ii)  uniform  distribution 
of  the  air  flow  along  the  collector  channels;  (iii)  the  physical 
properties  of  the  flowing  air  was  independent  of  temperature  and 
pressure;  (iv)  only  steady  states  were  considered  and  the  energy 
stored  in  the  components  was  neglected;  (v)  all  convective  heat 
transfer  coefficients  between  the  collector  components  and 
flowing  air  were  equal  and  constant  through  out  the  solution 
procedure  (vi)  the  temperature  of  the  flowing  air  was  uniform  over 
the  collector  cross-sectional  area.  The  energy  balance  equations  for 
the  collector  components  were  written  as  follows: 


For  the  absorber : 
\(  dTs 
Ac 


+  —  -kAn-x-  =  S 

A:V  dzjz= o 

For  the  plate  hrpb{Tp  -  Tb)  +  g-  (-Ank 


Ut(T  p  —  Ta)  +  hcpf(T  p  —  T  f)  +  hrpb(T  p  —  Tb) 

S  (93) 


dTs 

" s~dz 


z=D 


=  hcbf(Tb-Tf)  +  Ub(Tb-Ta) 


(94) 


For  the  metal  slats  : 


z— D 


csfLc(Ts  -  T f)dz 


z= 0 


(95) 


For  the  airflow  :  hcpf(T p  -  T f)  +  hcbf(Tb  -  T f) 

z=D 

+  T  y  2 hcs fnLc(Ts  -  T f)dz  =  qu  (96) 

z=0 

where  Ac  is  the  collector  total  surface  area  (WLc)  and  A  is  the  slat 
cross-sectional  area  (5Lc).  N  is  the  number  of  metal  vanes  per  meter 
of  width.  The  determination  of  the  radiative  heat  transfer  was 
given  in  Duffle  and  Beckmann  [52].  The  convective  heat  transfer 
coefficients  were  computed  on  the  basis  of  empirical  correlations 
developed  by  Heaton  et  al.  [21  ]  and  Nusselt  [53]  and  given  by  Ong 
[7].  The  solution  of  the  above  equations  was  carried  out  using  an 
iterative  procedure.  It  was  concluded  that  the  air  collector  with  a 
box-type  absorber  can  easily  be  modeled  mathematically  through 
the  development  of  expressions  for  its  parameters.  The  maximum 
number  of  metal  vanes  and  depth  of  the  air  duct  in  a  solar  air 
collector  are  interesting  factors  for  collector  design. 

The  mathematical  model  proposed  by  Ammari  [51  ]  was  similar 
to  the  above  for  the  energy  balance  of  the  absorber  plate,  back 
(bottom)  plate,  slats  metal  (vanes)  and  air  flow.  For  the  cover 
Ammari  [51]  has  proposed  the  follow: 

hnc(T  p  —  Tg)  +  hrpg(T  p  —  Tg)  =  hw(Tg  —  Ta)  +  hrga(Tg  —  Ta)  (97) 

The  sky  temperature  is  defined  as  [54]:  Tsky  =  0.05527^  5 .  The 
flow  is  assumed  to  be  hydro  dynamically  developed  at  the  collector 
inlet. 
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Kabeel  and  Mecarik  [50]  have  investigated  a  longitudinal  fins 
collector  with  two  covers.  This  model  is  show  in  Fig.  6c.  By 
considering  the  energy  balance  for  each  part  of  the  system  in 
steady  state  heat  flow  were  written  as  follow: 

For  the  top  glass  :  ctglS  +  hrgAg2{Tg2  -Tg i)  +  hcgxa\(Ja\  ~  Tgi) 

-  hw(Tgi  -  Ta)  =  0  (98) 


For  the  air  between  the  two  glasses  :  hcg2a\(Jg2  -  Ta i) 

+  ftcgiai  (Tg-[  -  Tai)  =  0  (99) 

For  the  second  glass  cover  :  Tg^ag2S  +  hcg2 f(Tf  -Tg 2) 

+  hrgip(T p  —  Tg2 )  -  hrg-[g2(Tg2  —  Tgi)  —  hcg2a-[  (Tg2  —  Tal)  =  0 

(100) 

For  absorber  plate  :  rgi  rg2apS  -  Zhpf(Tp  -T  f) 

-  hrpg2 (T P  -  Tg2)  -  Ub(Tp  -  Tfl)  =  0  (101) 

Z  is  the  ratio  between  the  actual  absorber  area  and  the  absorber 
area  normal  to  the  solar  radiation  falling  upon  the  collector. 

2.7.  Two  pass  solar  air  heaters 

The  concept  of  a  two-pass  air  heater  was  introduced  by 
Satcunanathan  and  Deonarine  [55]  and  later  considered  by  Caouris 
et  al.  [56]  for  liquid  based  systems.  These  authors  conducted 
experiments  on  the  two-pass  design  with  air  flow  between  the  two 
glass  covers  and  under  the  absorber.  Wijeysundera  et  al.  [57] 
motivated  by  a  need  to  consider  in  greater  detail  both  analytically 
and  experimentally  the  concept  of  two-pass  design,  developed 
heat  transfer  models  for  two-pass  flow  arrangements.  In  view  of 
uncertainties  which  exist  in  the  heat  transfer  coefficients  for 
various  processes,  they  are  made  the  following  assumptions  to 
keep  the  analysis  simpler:  (i)  The  temperature,  of  the  covers  and 
absorber  plate  were  uniform,  (ii)  For  the  thermally  developing  flow 
in  a  heater  channel,  suitably  averaged  heat  transfer  coefficients 
available  in  the  literature  were  used,  (iii)  The  bulk-mean 
temperature  of  the  cooling  fluid  was  changed  only  in  the  flow 
direction,  (iv)  The  area  of  the  collector  was  large  compared  to  the 
thickness  so  that  end  losses  and  shading  can  be  neglected,  (v) 
Effects  due  to  air  leakage  and  temperature  drop  in  the  plenums 
were  negligible,  (vi)  The  absorber  plate  and  the  covers  were 
diffuse-grey  for  long  wave  radiation.  For  incoming  solar  radiation 
the  covers  behave  like  specular  reflectors,  (vii)  The  collector  was  in 
quasi-steady  state.  Fig.  7  shows  a  schematic  view  of  the  two-pass 
flow  arrangement  under  consideration.  The  air  flows  first  between 
the  covers  and  then  through  the  bottom  channel.  There  is  a 
stagnant  air  gap  between  the  absorber  plate  and  the  cover  2.  For 
the  various  elements  of  the  collector,  Wijeysundera  et  al.  [57] 
wrote  the  nodel  energy  equations  as 

For  cover  1  :  =  Ac\aS  +  Achrg2g\(Tg2  -  Tgl)  -  Ac(hgla 

+  hw)(Tgl  -Ta)  (102) 


Cover  2 


For  cover  2  :  Q/2  =  A^S  -  Achrglg2(Tg2  -  Tgl) 

+  7\c(/irpg2  —  hcpg2)(Tp  —  Tg2)  (103) 

For  the  absorber  plate  :  Q.pj 

=  ACXPS  -  Achrpb(Tp  —  Tb)  —  Ac(hrpg2  +  hcpg2)(T p  —  Tg2)  (104) 


For  the  bottom  plate  :  Q,bf=  Achrpb(T  P-Tb)  -  Achcb  f(Th 


For  the  air  flow  between  covers 

=  Ac\J2\ 


Q/i 

h 


(T2  -  Tj)  +  F'(Tj  -  Tfl)  +  F'^^(r2  —  Tfl) 

n\f 


~Ta) 

(105) 

(106) 


For  the  airflow  between  absorber  plate  and  bottom  plate 


Q/2  =  AV21 


(Tg2  ~  )  +  F'(Tgl  -  Tfl)  +  F'  (Tg2  -Tfi) 

lira?  f 


cg2f 


(107) 


with  F'  =  mC/(hcg2 f  +  hcgl  f)WL[  1  -  exp(-(/icg2 f  +  hcgl  f)WL/rhC)}. 

Q  is  the  total  heat  flux  and  the  fractions  A,  of  the  solar  radiation 
absorbed  in  the  various  elements  of  the  collector  are  calculated 
using  the  net  method  [58].  A  similar  procedure  has  also  been 
applied  by  Wijeysundera  etal.  [57]  on  the  two-pass  solar  air  heater 
with  two  covers.  Wijeysundera  et  al.  [57]  used  this  mathematical 
model  to  study  the  thermal  performance  of  the  two-pass  solar  air 
heater  design.  The  results  show  that  the  performance  of  the  two- 
pass  design  with  air  flow  above  and  below  the  absorber  plate  in 
turn  is  better  than  the  single  pass  system.  The  results  of  the 
experiments  and  predictions  of  this  heat  transfer  model  have 
shown  that  there  is  some  scatter  in  the  experimental  data.  There 
results  were  also  compared  to  the  experimental  data  of 
Satcunanathan  and  Deonarine  [55]. 

A  new  mathematical  model  has  been  proposed  by  Verma  et  al. 
[59].  They  wrote  down  the  energy  balance  equations  for  the  design 
of  Fig.  8  under  usual  assumptions  [13]  to  show  the  optimization 
procedure  used.  The  two-pass  solar  air  heater  design  studied  by 
Verma  et  al.  [59]  is  different  from  the  one  modeled  by 
Wijeysundera  et  al.  [57]  (see  Figs.  7  and  8).  These  analyses  give 
the  following  equations: 


At  cover  glazing  :  agS  =  hrgs(Tg  -  Ts)  +  hcgn(Tg  -  T n) 

+  hw(Tg  —  T a)  +  hrgp(Tg  —  Tp)  (108) 


At  a  fluid  stream  between  cover  and  absorber  :  hcgf- 1 (Tg  -T^) 

=  Wr/,-rP)+^^f  d09) 

At  absorber  plate  :  aprgS  +  hrgp(Tg  —  Tp)  +  hCfip(T  fi  —  Tp) 

=  hcpf2(T p  -T f2)  +  hrpb(T p  -  Tb)  (HO) 


Fig.  7.  Schematic  view  studied  by  Wijeysundera  et  al.  [57]. 


Fig.  8.  Schematic  view  studied  by  Verma  et  al.  [59]. 
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At  a  fluid  between  absorber  plate  and  a  rear  plate  : 

^cp/2  (Tp  —  Tf2 )  =  hcf2b(Tf2  —  Tb)  +  f  —J— 


(111) 


At  a  fluid  stream  between  cover  and  absorber  :  hcgl(Tg  -  T/i) 

(120) 


.  _  ^  .  mC  f  dT  n  njr  ^  dT  n 

-ftc/ipOV,  -Tp)  =  — — +  M/C/-^- 


At  rear  plate  :  hcf2b(T f2  -  Tb)  +  hrpb{T P-Tb)  =  Ub(Tb  -  Ta) 

(112) 

The  differential  equations  can  be  solved  by  using  the  following 
boundary  conditions: 


r/1  Lorii  and  T/1  \X=L  ~  Tf2 


x=L 


(113) 


At  absorber  plate  :  apr gS  +  hrgp(Tg  —  Tp)  +  /ic/ip(^/i  —  Tp) 

-  hcpf2(TP  -  Tf2)  -  hrpb(Tp  -  Tb)  =  MpCp^f-  (121) 

At  a  fluid  between  absorber  plate  and  a  rear  plate  : 

hcpf2(Jp  ~  Tf2)  -hcf2b{Tf2  -  Tb)  =  +  M fC f 

(122) 


The  temperature  of  the  sky  Ts  have  been  calculated  by  using  the 
expressions  as  given  Swinbank  [54].  An  iterative  procedure  has  been 
used  by  Verma  etal.  [59]  to  solve  the  above  Eqs.  (109)-(113).  In  this 
procedure,  the  temperatures  of  the  cover,  plate,  air,  etc.  were  initially 
assumed.  The  new  temperatures  were  obtained  by  using  the 
solution  equations  and  were  compared.  The  procedure  was  repeated 
until  the  solution  is  converged  and  finally  the  collector  outlet  air 
temperature  was  obtained.  It  was  shown  that  there  exists  an 
optimum  mass  flow  rate  corresponding  to  an  optimum  flow  depth. 

The  mathematical  model  proposed  by  Wijeysundera  et  al.  [57] 
was  used  by  Naphon  and  Kongtragool  [17]  to  study  numerically  the 
effect  of  air  mass  flow  rate  on  the  heat  transfer  characteristic  and 
performance  of  the  single  two-pass  solar  air  heater.  The  explicit 
method  of  finite  difference  scheme  was  used  to  solve  these  models. 

In  order  to  study  the  effects  of  various  design  and  operational 
parameters  (duct  length  L,  depths  of  inlet  Di,  and  outlet  D2  air  flow 
channels  and  air  mass  flow  rate)  of  the  system,  Choudhury  et  al.  [33] 
modified  a  mathematical  model  proposed  by  Wijeysundera  et  al. 
[57].  The  schematic  view  is  identical  to  the  one  proposed  by  Verma 
et  al.  [59]  (see  Fig.  8).  The  steady-state  energy  balance  equations  for 
different  plates  and  flowing  air  in  the  inlet  and  outlet  channels  of  the 
air  heaters,  assuming  the  air  inlet  temperature  to  be  the  same  as  the 
ambient  temperature,  were  written  as  follows  [33]: 

For  the  two-pass,  single-cover  air  heater: 

Cover  glazing  :  agS  =  hw(Tg  -  Ta)  +  hcgfA  (Tg  -Tfl) 


+  hrgP{Tg  —  Tp)  (114) 

Fluid  between  cover  and  absorber  plate  :  mC(T f0  -  Ta) 

—  hcgf\(Tg  —  T  ft)  +  hcpfi(Tp  —  T  fj)  (115) 

Absorber  plate  :  aprgS  =  hcpfl(Tp  -Tfl)+  hcpf2(Tp  -  T f2) 

+  hrpb(T  p  —  Tb)  +  hrpg(Tp  —  Tg)  (116) 

At  a  fluid  between  absorber  plate  and  a  rear  plate  : 

mC(T0  -  Tf0)  =  hcgf2(Tp  -  T t2)  +  hcbf2(Tb  -  T p)  (117) 

At  rear  plate  :  hrpb(T P-Tb)~  hcb f2(Tp-T f2) 

=  Ub(Tb  -  Ta)  (118) 


At  rear  plate  hcf2b(T f2  -  Tb)  +  hrpb(T p  -  Tb)  -  Ub(Tb  -  Ta) 

=  MbCb ^  (123) 

The  initial  and  boundary  conditions  were:  T/i(x,0)  =  Ta(l), 
Tj2(x,  0)  =  Ta{  1),  Tp(x,0)  =  Tb(x,0)  =  Tg(x,0)  =  Ta{  1),  1^(0, t)  =  Ta  =  Tfil 
T/i(L,t)  =  Tf2{0,t),  Tf2(L,t)  =  Tf0. 

Choudhury  at  al.  [33]  have  investigated  a  double  pass  solar  air 
heater  with  double  cover.  This  model  is  shown  in  Fig.  9.  It  consists 
of  double  covers,  double  channel  design  with  double  air  flows 
between  inner  cover  and  absorber  plate  and  between  absorber  and 
bottom  plates  and  with  insulation  provided.  The  mathematical 


model  proposed  by  these  authors  is  as  follows: 

Cover  2  :  ag2S  =  hw(Tg 2  -  Ta)  +  /icg2gi  (Tg2  -  Tgl) 

+  drg2g\  (Tg 2  —  Tgi )  (124) 

Cover  1  :  ag\  rg2S  =  drg2g\  (T g\  Tg2)  +  hcg2g-i  (Tgl  Tg2) 

+  d-cg2fi  (Tgi  —  T f- 1 )  +  hrpg\  (Tgi  —  T p)  (125) 

Fluid  between  absorber  plate  and  a  cover  2  :  mC(T fl0  -  Ta) 

=  frcgi /i (Tgi  -  T/i )  +  fiCp/i  (Tp  —  T f\)  (126) 

Absorber  plate  :  oipxg2xg\S  =  hrpgl  (T p  -  Tgl)  +  hcpn  (Tp-Tn) 

+  dCp  f2  (T  p  —  T  /2)  +  Kpb(T  p  —  Tb)  (127) 

At  a  fluid  between  absorber  plate  and  a  rear  plate  : 

mC(Tf2o  -  T/io)  =  hcgf2(TP  —  7/2 )  +  hCbf2(Tb  ~  7/2 )  (128) 

At  rear  plate  hrpb(Tf2o  -Tb)~  hcbf2(Tb  -Tf2) 

=  Ub(Tb-Ta )  (129) 


Choudhury  et  al.  [33]  has  not  taken  account  of  the  heat  transfer 
by  the  radiation  from  the  cover  to  the  sky.  In  the  above  equations 
7}i  and  Tu  are  given  by:  Tn  =  (T/1o  +  Tfl)/2  and  Tfz  =  (T^o  +  T/o)/2.  The 
various  heat  transfer  coefficients  used  in  this  analysis  have  been 
computed  by  using  the  equations  given  in  Duffie  and  Beckann  [5], 
Tan  and  Charters  [6]  and  Hollands  et  al.  [61]  (see  Table  1). 


From  the  results  discussed  by  these  authors,  it  was  concluded 
that  the  optimum  values  of  the  duct  depths,  which  correspond  to 
the  minimum  annual  cost  per  unit  solar  energy  gain,  were  different 
for  different  duct  lengths  and  air  mass  flow  rates. 

Chauhan  et  al.  [60]  studying  the  drying  characteristics  of 
coriander  in  a  stationary  0.5  tone/batch  capacity  deep-bed  dryer 
coupled  to  a  double-pass  solar  single  air  heater  and  a  rock  bed 
storage  unit,  proposed  a  new  mathematical  model  for  this  air 
heater  design  (Fig.  8).  The  energy  balance  equations  for  different 
components  of  the  solar  air  heater  were  given  by  the  following: 

dT 

At  cover  glazing  :  MgCg-^  =  agS  -  hcgf]  (Tg  -  T/,) 

-hw(Tg-Ta)-hrgp(Tg-Tp)  (119) 


Fig.  9.  Schematic  view  of  energy-flow  of  double-pass  type  solar  air  heater  (see 
Choudhury  et  al.  [33]). 
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Recently,  to  evaluate  the  effects  of  two  kinds  of  roof  solar 
collector,  Zhai  et  al.  [18]  used  the  mathematical  model  developed 
by  Choudhury  et  al.  [33].  In  this  work,  for  the  natural  ventilation 
mode,  the  mean  air  temperatures  in  each  channel  were  calculated 
according  to  Hirunlabh  et  al.  [62]. 


T  fm  —  0.75  T20  +  0.25  Tr 


(130) 


and  the  air  mass  flow  rate  per  unit  of  the  system  was  given  by 


mfi  =  P 


An 
71  pi 


2gLsm(/3)(Ti0-Tr) 
Tr  ($in  +  Ho  +  f  b) 


(131) 


with  £,n  =  1.5,  £0  =  1.0 ,/=  0.056.  The  results  show  that  the  double 
pass  roof  solar  collector  is  more  efficient. 

Recently  Yamali  and  Solmus  [63]  have  investigated  a  desalina¬ 
tion  system  which  has  a  double-pass  flat  plate  solar  air  heater  with 
two  covers  as  a  component.  In  their  work,  a  transient  analysis  of 
double-pass  solar  air  heater  has  been  carried  out.  The  energy 
balance  equations  are  follows: 

dT 

For  cover  2  :  MgCg  =  agSWL  +  WL  hrgAg2  (Tg  1  ~Tg  2) 

—  WLhw(Tg2  —  Ta)  —  WLhrg2sky(Tg2  —  Ts/(y)  +  WL  hcgjg2(Tgj  —  Tg  2) 

(132a) 


with 

Tsky  =  Ta-6  (132b) 


Jt 

For  cover  1  :  MgCg  =  agzgSWLc  +  WL  hrpgX  (Tp  -  Tgl ) 

—  WL hcgi  (Tgi  —  Tfi)  —  WLhcgjg2(Tgi  —  Tg 2) 

-WLhcg,g2(Tg,-Tg2)  (133) 

dT 

First  air  pass  :  M fC f  =  WLhcpf-i  (Tp  -  T ^ ) 

+  WLhcgi  (Tg  1  -  T fi)  -  rrifCf(T f-io  -  T fU)  (134) 

Jt 

Absorber  plate  :  MgCg-^-  =  agr2gSWLc  -  WLhcpf2(Tp  -  Tf2) 

-  WLhcpf]  (T p  -  T„)  -  WLhrpg\  (T p  -  Tgl)  -  WLhrpb(Tp  -  Tb) 

(135) 


dT 

Second  air  pass  :  M fC f-gj-  =  WLhcpf2(Tp  -  T f2) 

+  WLhcb f2(Tb  -T f2)  —  m fC f(T f2o  -  T/l0)  (136) 

The  four  order  Runge-Kutta  method  was  used  to  solve 
numerically  the  energy  balance  equations.  For  this  purpose,  a 
computer  simulation  program  based  on  mathematical  model  was 
developed  by  means  of  MATLAB  software.  Theoretical  results 
illustrated  that  the  system  productivity  is  strongly  affected  by  the 
solar  air  heater  area  and  slightly  influenced  by  the  wind  speed 
variations  and  bottom  heat  loss  coefficients  of  the  solar  air  heater. 

With  the  aim  to  reduce  losses  from  the  front  cover  of  the 
collector,  Mohamad  [15]  modified  the  two-pass  solar  air  heater 
with  two  cover  proposed  by  Choudhury  et  al.  [33]  and 
Wijeysundera  et  al.  [57].  The  schematic  view  is  shown  in 
Fig.  10.  This  model  consists  of  double  covers,  double  channel 
design  with  double  air  flows  between  upper  cover  and  inner  cover 
and  between  inner  cover  and  absorber  plate,  the  lower  surface  of 
the  absorber  is  in  contact  of  the  insulation  provided.  Under  steady 
state  operating  conditions,  he  has  written  the  energy  balances  for 
this  design  as  follows: 

Cover  2  :  ag2S  =  hw(Tg2  -  Ta) +  hcg2fl(Tgi  -Tn) 

+  hrgig2  (Tg2  —  Tgi )  (137) 


Fluid  between  the  covers  :  hcg^  (Tg2  -  T f  t)  +  (T 

mC f  dT 
W  dx 


Tgi) 

(138) 


Cover  1  :  otg-[tg2S  —  hrg2g-[(Tg-[  Tg2)  +  hCgi  ji  (Tg\  Ty  i) 

+  A/2gi(Tgi  -  T f2)  +  hrgip(Tgi  -Tp)  (139) 


Fluid  between  the  cover  2  and  the  absorber  plate  : 

Tn  T*  dT 

hcf Tgl  (Tg,  -  r/2)  +  hcf2p(Tp  -  Tfl)  =  -£■  (140) 


For  the  absorber  plate  :  apxg2xg\S 

—  hc/2p(Tp  —  T  f2)  +  hrpg  i  (T  p  —  Tgi)  +  Ub(T  p  —  Ta) 


(141) 


In  order  to  increase  heat  transfer  from  the  absorber  to  the  air, 
Mohamad  [15]  suggested  replacing  the  absorber  plate  by  a 
porous  matrix  absorber.  In  this  case,  if  it  was  assumed  that  the 
air  and  the  solid  matrix  were  in  thermal  equilibrium  (the 
temperature  of  solid  was  equal  to  the  temperature  of  the  air 
locally),  the  energy  balance  in  the  second  passage  can  be 
simplified  as 


d2T 

h  /2gl  (Tgl  -  T  f2)  +  Ua  (T  a  -  T  f2 )  +  kef  f  dJ2  +  Ot  p  Tgl  rg2  S 

mC f  dT f2 
W  dx 


(142) 


The  results  indicate  that,  under  normal  operating  conditions, 
the  thermal  efficiency  of  the  suggested  heater  is  much  higher  than 
the  efficiency  of  the  conventional  air  heaters. 

2.8.  Double-pass  flat-plate  solar  air  heater  with  recycle 

Recently,  it  was  reported  that  increasing  the  fluid  velocity  by 
using  recycle  in  double-pass  parallel-plate  heat  exchangers 
enhances  the  heat-transfer  coefficient,  resulting  in  improved 
performance  [65].  The  schematic  view  of  the  design  of  the  solar 
heater  is  shown  in  Fig.  11.  In  addition  to  increasing  the  fluid 
velocity,  the  recycle  operation  also  produces  an  effect  by 
remixing  the  inlet  fluid  with  the  hot  outgoing  fluid.  Due  to 
the  fact  that  the  strength  of  forced  heat  convection  by  using 
recycle  effects  plays  an  important  role  in  designing  and 
operating  the  solar  air  heaters,  Flo  et  al.  [66]  developed  the 
mathematical  formulation  for  a  double-pass  solar  air  heater  with 
recycle.  In  this  model,  the  following  assumptions  were  made:  the 
temperatures  of  the  absorbing  plate,  bottom  plate  and  bulk  fluid 
are  functions  of  the  flow  direction  only.  Both  glass  covers  and 
fluids  do  not  absorb  radiant  energy.  The  radiant  energy  absorbed 
by  the  outlet  cover  is  assumed  to  be  negligible  [66].  With  the 
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heat  transfer  coefficient  is  constant  along  the  length  of  solar  air 
heater,  (iv)  thermal  conductivity  of  porous  media  is  constant  along 
the  length  of  solar  air  heater. 

The  mathematical  model  proposed  in  this  study  for  the  design  of 
the  double-pass  single  solar  air  collector  with  porous  media  is  based 
on  the  one  proposed  by  Naphon  and  Kongtragool  [17]  and 
Wijeysundera  et  al.  [57].  The  equations  of  the  top  glass  cover,  the 
first-pass  air,  the  absorber  plate  and  the  bottom  plate  are  similar  to 
given  by  Eqs.  ( 1 08),  ( 1 1 0)  and  (112).  With  the  presence  of  the  porous 
media  in  the  second-pass  air,  heat  transfer  Eq.  (Ill)  becomes: 


Fig.  11.  Schematic  view  of  energy-flow  of  double-pass  type  solar  air  heater  with 
recycle  (see  Ho  et  al.  [66]). 


above  assumptions,  the  energy  equations  were  written  as 
follows: 


Cover  1  :  hrpg2(Tp-Tgi)  +  h\{Th(z)-Tgl)  =  UgU{Tgi  -Ts) 

(143) 


Absorbing  plate  :  aprg i  rg2S  -  hT(J p  -  Ts)  +  hB(Tp  -  Ts) 

+  h,  ( Tp  -  Tb{x))  +  h2(Jp  -  Ta(x ))  (144) 

Bottom  plate  :  hrpR(J p  -  TR)  +  h'2(Ta(x)  -  TR) 

=  Ums(Tp-Ts)  (145) 


Channel  a  :  h2(T 
_  ~(R  +  1  )mCp 
W 


p  -  T Q(X))  -  h'2(T a(x) 
dTq(x) 
dx 


(146) 


Channel  b  :  h^(Tp  -  Tb(x))  -  h\  (Tb(x)  -  Tgl) 


{R  +  1  )mCp 
W 


dTh(x ) 
dx 


(147) 


In  this  equation,  the  characteristic  length  is  the  equivalent 
diameter  of  the  duct.  For  the  open  conduit,  lower  channel  and 
upper  channel,  they  are  respectively  given  by: 


4HW  n  _  (4AH)W 

e’°  -  2 H  +  W ’  De  a  ~  (2  AH)  +  W  and 

[4(1 -A  )H]W 
e,a  ~  [2(1  -  A)H]  +  W 


(148) 


and  the  air  velocities  of  the  open  conduit,  lower  channel  and  upper 
channel  used  by  Ho  et  al.  [66]  are  given  by 


m  _  (1  +R)m  _  [m(l+/?)] 

=  (2 HW)p’  Va  =  (2A H)Wp  3n  Vb  =  [2(1  -  A )H]Wp 

(149) 


A  is  the  ration  of  channel  thickness. 

The  experiment  had  been  carried  out  by  these  authors.  With  an 
iterative  procedure,  Ho  et  al.  [66]  investigated  the  effect  on  the 
collector  efficiency  with  ratio  as  the  variable  parameter  for  this 
type  of  solar  air  heater  design  (see  Fig.  11).  The  results  showed  that 
the  theoretical  prediction  agree  reasonably  well  with  experimental 
data.  It  was  also  confirmed  that  the  collector  modules  with 
external  recycle  are  thermo-hydraulically  better.  This  is  economic¬ 
ally  feasible  in  the  design  of  double-pass  flat  solar  air  heaters. 

2.9.  Double-pass  flat  plate  solar  air  heater  with  porous  media 

Naphon  [69]  studied  theoretically  the  heat  transfer  character¬ 
istics  and  performance  of  double-pass  flat  solar  air  heater  with  and 
without  porous  media.  The  model  was  based  on  that  of  Naphon 
and  Kongtragool  [17]  with  the  following  assumption:  (i)  flow  of  air 
is  steady,  (ii)  outside  convective  heat  transfer  coefficient  is 
constant  along  the  length  of  solar  air  heater,  (iii)  inside  convective 


kp  is  the  thermal  conductivity  of  the  porous  media.  The  implicit 
method  of  finite-difference  scheme  had  been  employed  to  solve 
these  models.  The  effect  of  the  thermal  conductivity  of  the  porous 
media  on  the  heat  transfer  characteristics  and  performance  have 
been  considered.  The  theoretical  results  obtained  from  the  model 
were  validated  by  comparison  with  the  experimental  data  of 
Sopain  et  al.  [70]  and  there  was  reasonable  agreement. 

2.10.  Double-pass  flat  plate  solar  air  heater  with  longitudinal  fins 


Naphon  [71]  has  investigated  heat  transfer  characteristics  and 
entropy  generation  of  the  double-pass  flat  plate  solar  air  heater 
with  longitudinal  fins.  The  basic  physical  equations  used  to 
describe  the  heat  transfer  characteristics  of  the  solar  air  heater,  as 
shown  in  Fig.  12a  and  b,  are  developed  from  the  conservation 
equations  of  energy. 

With  the  first  law  of  thermodynamics  and  based  on  the  Naphon 
and  Kongtragood’s  model  [17],  Naphon  developed  a  mathematical 
model  for  double-pass  flat  solar  air  heater  with  longitudinal  fins. 
Using  the  assumption  proposed  in  the  case  of  double  pass  solar  air 
heater  with  porous  media,  the  following  equations  have  been 
developed: 

For  the  glass  cover  :  agS  =  hw(Tg  -  Ta)  +  hCf-[Cg(Tg  T /i) 

+  hrgp(Tg  —  T  p)  +  hrag(Tg  —  T  a)  (151) 


dT 


For  the  first-pass  air  :  mC =  hCf\g(Tg  -T f- 1) 

N  rz=H 

+  hCfiP(T p  —  T fi )  +  -  /  2Lh\  (7V \  —  T f\)dz 

Afrontal  Jz= 0 


(152) 


Fig.  12.  (a)  Schematic  diagram  of  the  solar  heater  (Naphon  [71]  model),  (b) 
Schematic  diagram  of  fin  section  (Naphon  [71]  model). 
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Fig.  13.  (a)  Three  pass  single  cover  air  heater  (Chandra  et  al.  [72]).  (b)  Three  pass  two  covers  air  heater  (Choudhury  et  al.  [75]). 


where  m  is  the  mass  flow  rate  of  fluid  per  unit  with,  j4frontai  is  the 
frontal  area  of  the  channel,  H  is  the  height  of  the  fin,  z  is  the 
distance  from  the  absorber  plate,  Tvl  is  the  fin  temperature  in  the 
upper  channel  and  L  is  the  length  of  the  collector.  N  is  the  number 
of  fins. 


For  the  upper  fin  : 


N 


'Z=H 


A 


N 


^frontal 


-kA 


dT 


v\ 


dz 


frontal  Jz= 0 


z= 0 


2Lhi  (Tv i  —Tfi  )dz 


(153) 


For  absorber  plate  :  agx gS  =  hcfip(T ,,  -  T  f1)  +  hcf2p{T p  -  T f2) 


+  hrgp(Tp  —  Tg)  +  hrpb(Tp  -  Tb)  +  — - 

-'‘frontal 


dT\n\ 
dz  ) 


z= 0 


+ 


dTv2\ 
dz  ) 


(154) 


2.2  2.  Triple  pass  solar  air  heaters 

The  schematic  configuration  of  a  typical  triple-pass  solar  air 
heater  is  shown  in  Fig.  13a.  It  consists  of  one  cover,  an  absorber 
plate  and  two  rear  plates  (Fig.  13a).  Chandra  et  al.  [72]  have 
investigated  this  design.  Neglecting  the  thermal  capacity  of  the 
plates,  the  energy  balance  equations  were  written  as  follows: 


Cover:  agS  =  hrgs(Tg  -  Ts)  +  hw(Tg 
+  frrgpl  (Tg  —  T /3) 


Ta)  +  hcgf3(Tg  -  T f3) 

(158) 


TT1  Hi ' 

FluidB:  -^l-Jl  =  hcgf3(Tg-Tf3)-hCplf3(Tf3-TpX)  (159) 


Tv 2  is  the  fin  temperature  in  lower  channel  and  Asf  is  the  cross- 
section  area  of  fin. 

-z=H 

I  V  I 

For  lower  fin  : 


N 


A 


-kA 


frontal 


-^frontal 

dTy  2 
dz 


[  2Lh2(Tv2  -  T f2)dz 

Jz= o 


z= 0 


(155) 


For  second-pass  air :  mC 


dTf2 


f  dx 

=  hCf2P(T p  -  T f2)  +  hCf2p(Tp  -  T f2) 

N 

+ 


-^frontal  Jz= 


[Z  H  2Lh2  (Tv2  —  T f2)dz 
Jz=0 


(156) 


For  bottom  plate  :  hcf2b(Jb  -  T f2)  +  hrpb(T p  -  Tb) 

+  Ua(Tb-Ta)  =  0  (157) 

It  was  found  that  the  thermal  efficiency  increases  with 
increasing  the  height  and  number  of  fins. 


Plate  1  :  TgCXpi S  =  hrgv\  (Tg  -  Tpi)  +  /icpl /3(Tpi  -  T /3) 

+  hrp\r2(T pi  -  T p2)  +  hcpif2(Tp\  -  T f2)  (160) 

_ ^  HT 

Fluid  2  :  —  f  =  hcpl  f2(T p\  -  T f2)  -  hcp2f2(T f2  -  T p2) 


(161) 

Plate  2 

:  ^rpi p2 (T pi  -  T p2)  +  hcf2p2(T f2  —  Tp2) 

=  hrp2b(Tp2  -  Tb)  +  hcp2fi(Tp2  —  Tfi) 

(162) 

Fluid  1 

■■  nwLdir  =  hcPfi  ( Tp2  -  Tfi)  -  hcbf ,  (Tfi  -  Tb) 

(163) 

Plate  3 

-  Kp2b(T p2  T p3)  +  hcbfi  (Tfi  -  Tb)  =  Ur(Tb  -  Ta ) 

(164) 

with  7/i (x  =  0)  =  Ta,  Tn(x  =  L)  =  TP(L)  and  Tp( 0)  =  T&( 0). 

In  order  to  appreciate  the  analytical  results,  numerical 
calculations  have  been  carried  out  corresponding  to  the  experi- 
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mental  data  of  Ezeike  [74].  The  results  of  this  theoretical  work 
differ  substantially  from  Ezeike’s  results.  The  results  obtained  from 
an  exact  analysis  and  from  an  approximate  equivalent  thermal 
network  method  were  almost  consistent. 

This  model  of  solar  air  heater  has  also  been  investigated  by 
Choudhury  et  al.  [75].  They  assumed  that  the  air  inlet  temperature 
is  the  same  as  the  ambient  temperature.  In  the  mathematical 
model  proposed  by  these  authors,  Eq.  (158)  is  reconsidered 
without  the  term  of  radiative  exchange  between  the  cover  and  the 
sky;  Eqs.  (160),  (162)  and  (163)  are  the  same.  Eqs.  (159),  (161)  and 
(163)  have  been  replaced  by  the  following: 


Fluid  3  :  mCf{T f0-T f20) 


—  hcgf3(Tg  T  f3)  hcpif3(Tf3  Tpi) 

(165) 

Fluid  2:  mCf(Jf2o-Tno) 

=  ^cpl  /2  (7pl  -  T  f2)  -  hcp2f2{T  f2  -  7p2) 

(166) 

Fluid  1  :  mCf(Tfl0-Ta) 

=  frcp2/i  {Tp2-Tn)~  hcp3f-[  (T/!  -  T  p3) 

(167) 

Another  three-pass  solar  air  heater  was  investigated  by 
Choudhury  et  al.  [75].  It  consists  of  two  covers,  an  absorber  and 
a  rear  plate  (Fig.  13b).  Assuming  the  air  inlet  temperature  to  be  the 
same  as  the  ambient  temperature,  the  steady  state  energy  balance 
equations  for  the  different  plates  and  flowing  air  in  the  inlet, 
intermediate  and  outlet  channels  of  the  air  heaters  can  be  written 


as  follows: 

Cover  1  :  ag^S  =  hgia(Tgi  -  Ta)  +  ^rg2gl  (Tgl  Tg2) 

+  frcgi /i  (Tgi  -  T fi )  (168) 

Fluid  1  :  mCf(T fl0  -  Ta) 

=  frcgi/i (Jg\  -  T fi)  -  hcg2f\ {Tg 2  -  T fi)  (169) 

Cover 2  :  xglap2S  =  /icg2/1(Tg2  —  T n)  +  /icg2/2(Tg2  -  T f2) 

+  frrg2gl  (Jg2  —  Tgi)  +  hrpip2(Tpi  —Tp2)  (170) 

Fluid  2  :  mC(J f2o  -  T no) 

=  hcg2f2{Tg2  —  Tf2)  —  hcpi/2(^pi  -  T f2)  (171) 


Absorber  :  xg\Xg2ap\S  =  hcpX /2(Tpi  -  T f2)  +  hcpl f3(T P\-T f3) 


+  ^rpig2 (7pi  -  Tg2)  +  hrplp2(Tp1  -  Tp2) 

(172) 

Fluid  3  :  mCf(Tf0-Ta) 

=  ^cpi /3 (T pi  -  T f3)  —  hcp2f3(T f3  —  T p2) 

(173) 

Rear  plate  :  hrpip2(T pi  -  T p2)  +  hcp2f3(T f3  -Tp2) 

=  Ur(Tp2  -  Ta) 

(174) 

In  all  these  equations  Tf i,  Tp  and  Tp  are 

defined  as 

l/i  =  (Tfio  +  Ta)l 2,  Tf2  =  ( Tj2o  +  T/i0)/2  and  Tf 3  =  ( Tfo  +  TJf2o)/2,  respec¬ 
tively.  In  this  study,  the  authors  have  made  an  effort  to  find  the 
effects  of  duct  length,  duct  depths  of  inlet  (Di),  intermediate  (D2) 
and  outlet  (D3)  air  channel  and  mass  flow  rate  (m)  on  the  cost 
benefit  ratio.  Efforts  were  also  made  to  identify  the  most  superior 
air  heater  by  comparing  the  cost  effective  performance  of  these  air 
heaters  with  those  of  single  pass  with  no,  single  and  double  cover 
and  two  pass  with  single  and  double  covers  air  heaters  studied 
earlier  by  the  authors  [33,75].  Choudhury  et  al.  [75]  in  their 
analysis  have  used  the  heat  transfer  coefficients  predicted  by  the 
equations  given  in  Duffie  and  Beckmann  [5],  Tan  and  Charter  [6] 
and  Hollands  and  Shewen  [61  ].  It  was  shown  that  the  double  cover, 
three  pass  solar  air  heater  with  air  flow  from  top  to  bottom  is 
observed  to  be  more  cost-effective  than  the  single  cover,  three  pass 
solar  air  heater  with  air  flow  from  bottom  to  top. 

2.22.  Multi-pass  flat-plate  solar  air  heaters  with  external  recycle 

Ho  et  al.  [76]  have  developed  a  theoretical  formulation  for  a 
multi-pass  solar  air  heater  with  external  recycle  and  investigate 
the  recycle  effect  on  collector  efficiency.  The  effects  of  channel 
thickness  ratio  on  the  collector  efficiency  was  also  studied.  In  this 
study,  Ho  et  al.  [76]  have  considered  a  recycled  four-pass  solar  air 
heater  with  vertical  thickness  2aH  and  2(1  -  a)H  (a  is  the  channel 
thickness  ratio  in  vertical  direction),  and  horizontal  width  2/3H  and 
2(1  -  fl)H  (/3  is  the  channel  thickness  ratio  in  horizontal  direction), 
respectively,  which  was  obtained  by  inserting  an  absorbing  plate 
and  insulation  sheet  with  negligible  thickness  into  a  parallel 


Fig.  14.  (a)  Multi-pass  flat-plate  solar  air  heater  with  recycle  (Ho  et  al.  [76]).  (b)  Schematic  view  of  energy-flow  of  multi-pass  type  solar  air  heater  with  recycle  (see  Ho  et  al. 
[76]). 
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conduit  of  height  2H,  length  L,  and  width  W  (see  Fig.  14).  Assuming 
that:  (i)  the  temperatures  of  the  absorbing  plate,  bottom  plate  and 
bulk  fluid  were  functions  of  the  flow  direction  only,  (ii)  both  glass 
covers  and  fluids  do  not  absorb  radiant  energy,  (iii)  the 
absorptance  of  the  glass  cover  2  was  negligible,  the  energy 
equations  have  been  written  as 

Cover  1  :  hrpgi  (Tp  -  Tgl )  +  h'3(Tc(x)  -  Tgl)  +  h'4(Td(x)  -  Tgl ) 

=  Ugls(Tgl  -  Ts)  (175) 


Absorbing  plate  :  apTgxT g2S  -  hT{Tp  -  Ts)  -  hB(T  p  -  Ts) 

-  hi  (Tp  -  Ta(x))  -  h2(Tp  -  Tb(x))  +  h3(Jc(x )  -Tp)  +  h4(Td(x) 

—  Tp)  =  0  (176) 


Bottom  plate  :  fi[(Ta(x)  -  TR)  +  h'2(Tb(x)  -  TR)+hrpR(Tp-  TR) 

=  UBs(Tp-Ts)  (177) 


Channel  a(flowing  under  the  absorber  plate) 

hi(Tp  —  Ta(x))  -  h\  (Ta(x)  -Tr) 


(R  +  l)mCp 

Wp 


dTa(x) 

dx 


(178) 


Channel  b  :  h2(Tp  -  Tb(x))  -  h2(Tb(x)  -  TR) 


-(R  +  l)mCp 


(1  -p)W  \  dx 


dTb(x ) 


(179) 


Channel  c  :  h3(Tp  -  Tc(x))  -  h'3(Tc(x)  -  Tgl) 

_  '-RmCp  1  dTb(x) 

=  Wp  \  ~dx~ 


(180) 


Channel  d  :  h4(Tp  -  Td(x))  -  h4(Td(x)  -  Tgl) 

_  '  RmCp  1  dTb(x) 

=  (1  -  P)W\  ~dx~ 


(181) 


In  this  work,  the  equivalent  diameters  were  defined  as 


(4aH)WP 

ep  (2  aH)  +  WP  ’ 

[4(1  -a)H]WP 

e'c  [2(1  —  a)H]  +  WP’ 


[4aH\W(\  -  P) 

[2 aH]  +  W(  1  -  P)  ’ 

[4(1  -  o')H]W(l  -  p) 
ed  [2(1  —  <x)H\  +  W(1  —  P) 

(182) 


and  average  air  velocities  given  by  the  following: 


Vr 


m, 


Vr 


(2aH)W/3p 5 
mc 


Vb 


2H(1  -aW/Sp1 


mb 

2HaW(\  -p)p' 

-  _  md 

d~2H(1  -a)W(  1  -  P)p 


(183) 


The  above  set  of  equations  was  solved  numerically  using  an 
iterative  procedure.  In  this  procedure,  the  prescribed  collector 
geometries,  system  properties  and  operating  conditions  were 
given,  as  well  the  initial  guesses  of  the  mean  temperatures  were 
specified  for  calculating  the  physical  properties.  The  new  values  of 
mean  temperatures  were  then  re-calculated.  If  the  re-calculated 
average  temperature  values  were  not  quite  close  to  the  assumed 
values,  the  iteration  calculations  performed  above  will  be  done 
until  the  last  assumed  values  meet  the  final  calculated  values.  The 
results  have  shown  that,  the  introduction  of  recycle  effect  was  a 
feasible  way  to  increase  the  collector  coefficient  due  to  the 
increasing  fluid  velocity  effect  and  may  compensate  for  the 
temperature  difference  decrement. 


2.13.  Solar  air  heater  with  a  compound  parabolic  concentrator 

Fig.  15  shows  this  solar  air  model.  Pramuang  and  Exell  [77] 
reported  the  results  of  an  experimental  study  in  which  the  method 


Fig.  15.  Schematic  view  of  a  solar  air  heater  with  a  compound  parabolic 
concentrator  (Premuang  and  Exell  [77]  and  Tchinda  [78]). 


of  Chungpaibulpatana  and  Exell  is  used  to  determine  the  collector 
parameters  of  a  solar  flat  plate  collector  with  a  CPC  for  heating  air. 
Recently,  Tchinda  [78]  has  quantified  the  heat  transfer  within  this 
design  of  solar  air  heater  and  a  mathematical  model  analysing  the 
collector  thermal  performance  has  been  proposed.  In  this  model 
(see  Fig.  15),  the  following  assumptions  were  made:  (i)  It  was 
assumed  that  the  CPC  is  ideal  and  free  from  fabrication  errors,  (ii) 
Any  beam  of  radiation  incident  within  the  acceptance  angle  0a, 
with  the  help  of  the  parabolic  reflector  can  reach  the  receiver,  (iii) 
The  reflection  of  radiation  from  the  parabolic  reflector  was 
accounted  for  by  the  apparent  reflectance  p<n> .  (iv)  The  direction 
of  the  beam  radiation  incident  on  various  components  in  the 
collector  can  be  found  though  geometry,  (v)  The  succeeding 
absorption  and  transmission  processes  inside  CPC  were  diffusive 
and  were  accounted  for  in  terms  of  diffuse  properties.  The  solar  and 
infrared  energy  exchanges  in  the  collector  were  treated  separately 
using  pertinent  radiative  properties  in  the  spectrum,  (vi)  Physical 
and  optical  properties  of  materials  were  assumed  to  be  indepen¬ 
dent  of  temperature,  (vii)  The  concentrator  did  not  produce  an 
image  of  the  light  source;  hence  it  was  called  non-imaging 
concentrators.  With  the  above  assumptions,  the  energy  equations 
may  be  written  as  follows: 

A 

For  cover  :  McMpc-^  =  /(t)[ac  +  acicppPmn)]j-+  hRp(T p  -  Tc) 
+  hp/c(Tp  -  Tc)  -  hgsFc  -  Ts)  -  hc/a(Tc  -  T„)  (184) 


with  t  >  0. 


For  plate  absorber  :  MPC 


m?cp$p 

<J  u(9 


dTP 
p^pp~dt 
A, 


Dtp  +  OipppPc  ^ 


A 


hRp(Tp  —  Tc)  —  hp/c(T  p  —  Tc) 

(185) 


with  t  >  0. 


For  fluid:  pfefCpf^f 

=  vM.^.Uo(Tf.Tb)  (186) 

with  t  >  0  and  0  <x<L.  where  P  =  1  -  g/lp  ([80]),  g  is  the  gap 
thickness.  AC=W*L,  Ap  =  lp*L  and  qu  =  Uj(Tp  -  Tf).  An  iterative 
method  was  used  by  the  author  to  bring  the  effect  of  the 
temperature  dependence  of  various  heat  transfer  coefficients.  For 
certain  temperatures,  they  were  first  calculated  by  using  the 
standard  expressions  given  earlier.  Equations  were  solved  by 
assuming  the  constant  and  the  new  solutions  used  to  generate  all 
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Fig.  16.  (a)  v-Groove  air  collector,  Karim  and  Hawlader  [83 ].  (b)  v-Groove  air  collector  (Bashria  et  al.  [84]).  (c)  v-Groove  air  collector  with  porous  media  (Bashria  et  al.  [85]).  (d) 
v-Groove  air  collector  (Kabeel  and  Mecarik  [50]). 


the  heat  transfer  coefficients  again  till  the  values  converge.  This 
model  has  been  validated  with  the  experimental  data  of  Pramuang 
and  Exell  [77]. 

2A4.  V- groove  solar  air  collector 

The  design  given  in  Fig.  16a  was  also  studied  experimentally 
and  theoretically  by  Karim  and  Hawlader  [82,83].  The  energy 
balances  for  absorber  plate,  back  insulation  and  working  fluid  of 
the  above  collector  are  given  by  the  following  Karim  and  Hawlader 
[83]: 

rgapS  +  Ut(Ta  -  Tp)  +  hi  (T f  -  Tp)  +  hr(Tb  -Tp)  =  0  (187) 

Ub(Ta  -  Tb)  +  h2(Jf  -  Tb)  +  hr(Tb  -Tp)  =  0  (188) 

hi  (T  p  —  T  f)  +  h2(Tb  —  T  f)  =  qu  (189) 

A  series  of  experiments  were  conducted,  data  were  recorded  for 
different  operating  variables  to  determine  the  performance  of  v- 
groove  solar  air  collectors.  Experimental  and  analytical  results 
have  shown  a  good  thermal  performance  of  the  v-groove  collector. 
Satisfactory  qualitative  and  quantitative  agreement  between 
experimental  and  analytical  results  was  achieved. 

An  Internet  based  mathematical  simulation  has  been  conducted 
and  developed  by  Bashria  et  al.  [84,85]  and  Bashria  et  al.  [86]  to 
predict  the  thermal  performance  for  different  designs  of  solar  air 
heaters  with  v-groove  absorbers,  single,  double  or  triple  glass 
covers. 

Fig.  16a-c  show  the  designs  of  v-groove  solar  air  collector 
studied.  A  number  of  simplifying  assumptions  has  been  made  to 
lay  the  foundation  without  obscuring  the  basic  physical 
situation.  These  assumptions  are  as  follows:  (i)  performance  is 
steady  state;  (ii)  There  is  no  absorption  of  solar  energy  by  a  cover 
insofar  as  affects  losses  from  the  collector,  (iii)  Heat  transfer  fluid 
is  considered  a  non-participating  medium,  (iv)  The  radiation 


coefficient  between  the  two  air  duct  surfaces  is  found  by 
assuming  a  mean  radiant  temperature  equal  to  the  mean  fluid 
temperature,  (v)  Loss  through  front  and  back  are  of  the  same 
temperature.  The  steady  state  energy  equations  yield  the 
following  equations: 

In  the  single  pass  v-groove  absorber: 

For  the  cover  :  Ut(Tg  -  Ta)  =  hi(Tp  -  Tg)  +  hr(Tp  -  Tg)  (190) 

^-groove  absorber  :  rgapS  =  hi  (Tp  -  Tg)  +  hr(Tp  -  Tg) 

+  h2(T p  -  T f)  (191) 

Fluid  medium:  h2(Tp  -  T  f)  +  Ub(Ta  -  Tf)  =  ^  (192) 

In  double  pass  double  duct: 

For  the  cover  :  Ut(Tg  -  Ta)  =  hi(Tp  -  Tg)  +  hr(Tp  -  Tg)  (193) 

^-groove  absorber  :  rgapS  =  h2  (Tp  -  T/i )  +  hi  (Tp  -  Tg) 

+  hr(T  p  —  Tg)  +  h2(Tp  —  T  f2)  +  hr2(T  p  —  Tr) 

Fluid  medium  in  the  upper  passage  :  h2(Tp  -  T/i) 

mCAdTf i 
W  )  dx 


(194) 

(195) 


Fluid  medium  in  the  lower  duct :  li3(Tp  -  T f2)  +  h4(T S2 
( mC  f  \  dT  f2 
=  \~W~ )  ~dx~ 

Bottom  plate  :  Ub(Tr  -  Ta)  =  h4(T f2  -  Tr)  +  hr2(Tp  -  Tr) 


~Tr) 

(196) 

(197) 


In  double  flow  double  duct  packed: 


For  the  cover:  Ut(Tg  -Ta)  =  h,(Tp  -  Tg)  +hr(Tp  -  Tg)  (198) 
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^-groove  absorber  :  rgapS  =  h2  (Tp  -  T  fX )  +  h\  (Tp  -  Tg) 

+  hr(Tp  -  Tg)  +  h3(Tp  -  Tf2)  +  hr2(Tp  -  Tpr)  (199) 


Fluid  medium  in  the  upper  passage  :  h2(Tp  -T 

_  fmCpX\  dT fi 
(  w  )  dx 


(200) 


Fluid  medium  in  the  lower  duct : 


+  h5(r/2  -Tr) 


( mCp2\  dT f2 
ylAT )  ~dx~ 


h3  (T  p  —  T  f2)  +  h4(T  f2  —  Tpr) 

(201) 


Porous  media  :  h4(J pr-T f2)  =  hr2(T p  -  T pr)  +  h6(T pr  -  Tr) 

(202) 


Bottom  plate  :  U5{T f2  -  Tr)  +  h6(T pr  -  Tr)  =  Ub(Tr  -  Ta)  (203) 


The  pressure  drop  through  the  channel  in  the  air  heater  has 
been  computed  using  the  following  expression  [87]: 


P  = 


fo+yL 


(204) 


where  for  f0  =  24/Re,  y  =  0.9  for  laminar  flow  (Re  <  2550), 
/o  =  0.0094,  y  =  2.92 Re~° 15  for  transitional  flow  (2550  <  Re  < 
104)  and  f0  =  0.059  Re~02,  y  =  0.73  for  turbulent  flow 
(104  <  Re  <  105). 

This  mathematical  model  has  been  used  recently  by  Bashria 
et  al.  [86]  to  find  the  effect  of  different  parameters  such  as  mass 
flow  rate,  flow  channel  depth  and  collector  length  on  the  system 
thermal  performance  and  pressure  drop  through  the  collector  (see 
Fig.  16a-c).  The  prototype  internet  program  was  developed  to  be 
used  as  a  tool  to  support  the  design  of  solar  air  heaters  used  for 
drying.  The  results  were  validated  by  comparison  the  predicted 
output  and  the  experimental  results  carried  out  by  Karim  and 
Hawlader  [82].  A  great  correlation  has  occurred  between  the 
experimental  and  the  predicted  efficiencies.  Also  a  good  correla¬ 
tion  has  been  found  between  the  experimental  and  predicted 
efficiencies  and  outlet  temperature,  respectively. 

Kabeel  and  Mecarik  [50]  have  investigated  v-groove  air  heater 
system  with  two  covers  and  the  flow  passes  over  and  under  the 
absorber  plate  (Fig.  16d).  By  considering  the  energy  balance  for 
each  part  of  the  system  in  the  steady  state  heat  flow  were  written 
as  the  follow: 

For  the  upper  glass  cover  :  aglS  +  hrg\g2(Tg2  -  Tgl) 

+  frcgiai  (Ta  i  —  Tgi )  —  hw(Tgi  —  T  a)  —  0  (205) 


For  the  second  glass  cover  :  xg\  ag2S  -  hcg2a i  ( Tg2  -  Ta i) 

—  h2(T  f2  —  Tg2)  —  hrg  2gi  (Tg  2  —  Tgi)  +  hrg2p(Tp  —  Tg2 )  —  0 


(206) 


For  the  flow  between  the  two  glass  covers  : 

—  ftglal  (Ta l  —  Tgi )  =  0 


dg2a\  (Tg2 


Tal) 

(207) 


For  the  absorber  plate  :  rgi rg2apS  -  Zhp fX  (Tp  -Tfi) 

—  Zhpf2(T  p  —  T  f2)  —  hrpg2(T  p  —  T g2)  —  hrpb(T  p  —  Tb)  —  0  (208) 

For  Back  plate  :  hrpb(T P-Tb)  +  hcb n  (T n  -  Tb)  -  Uh(Th  -  Ta) 

=  0  (209) 

Z  is  the  ratio  between  the  actual  absorber  area  and  the  absorber 
area  normal  to  the  solar  radiation  following  upon  the  collector. 
With  this  mathematical  model  taking  into  account  the  absorber 
shape  factor,  Kabeel  and  Mecarik  [50]  have  outlined  the  effect  of 
the  change  of  the  triangular  angle  0  on  the  pressure  drop,  on  the 
heat  transfer  added  to  the  flowing  air  and  deduced  the  optimum 
angle  of  the  v-groove. 


2.25.  Single  pass  solar  air  heater  with  packed 
2.15A.  Type  2 

The  type  1  solar  heater  is  shown  in  Fig.  17a.  It  consists  of 
conventional  configuration  of  Fig.  4  with  the  flow  passage  filled 
with  packing.  In  the  mathematical  model  proposed  by  Choudhury 
and  Garg  [88],  energy  balance  equation  of  the  cover  is  given  by 
Eq.  (36).  For  the  absorber  plate,  packing  and  air  flow  and  bottom 
plate,  energy  balance  equations  in  the  steady  state  were  given  by 
the  follow: 

Absorber  plate  :  rgapS  =  (hrgp  +  hCgp)(Tp  -  Tg) 

+  drpp-[(T  p  —  T  pi)  +  hcpf(T  p  —  T  f)  (210) 

Packing  :  hrpp\  (Tp  -  Tpl )  =  hcpl f(T pi  -  T f)  +  hrplb(T pi  -  Tb) 

(211) 

TT\  ^  dT 

Air  flow/ :  =  hcpf(Tp-Tf)  +  hcpU(Tin  -Tf) 

+  hcbf(Tb-Tf)  (212) 


Bottom  plate  :  hrpW(T pX  -  Tb)  =  hcb f(Tb  -Tf)  +  Ub(Tb  -  Ta) 

(213) 

The  boundary  conditions  for  the  above  design  are: 
Tj(x  =  0)  =  Tfi,  Tj(x  =  L)  =  Tfo.  In  this  model,  it  was  assumed  that 
the  heat  transfer  resulting  from  conduction  in  the  packed  flow 
channels  is  neglected  due  to  point  contacts  between  packing  and 
the  plates. 

Recently,  Paul  and  Saini  [90]  have  used  the  mathematical  model 
proposed  by  Choudhury  and  Garg  [88]  to  outline  the  bed 
parameters  in  such  a  way  that  the  solar  energy  collection  systems 
deliver  energy  with  minimum  cost.  A  computer  program  in  ‘C 
language  has  been  developed  to  obtain  thermal  efficiency  and  cost 
per  unit  energy  delivered  by  packed  bed  solar  air  heater  having 
wire  mesh  screens  and  pebbles  as  packing  materials  with  different 
specifications.  This  program  was  validated  by  comparing  some  of 
the  predicted  values  of  thermal  efficiency  with  experimental 
values  obtained  by  Collier  [92]. 

2.25.2.  Type  2 

It  consists  of  two  sheets  of  cover  and  the  back  plate;  the 
passage  between  the  inner  cover  and  the  back  plate  is  provided 
with  black-painted  packing,  which  act  as  absorbers  of  solar 
radiation  and  transfer  the  heat  to  the  air  flowing  through  the 
packed  channel  (see  Fig.  17b).  This  device  has  been  investigated 
by  Choudhury  and  Garg  [88]  and  the  energy  balance  written. 
Energy  balance  of  the  outlet  cover,  back  plate  and  air  flow 
were  given  by  Eqs.  (36),  (212)  and  (213)  respectively.  For 
the  inner  cover  and  packing,  the  energy  balances  are  as 
follows: 

For  the  inner  cover  :  xg2ag\S  +  hrp^g2(Tp-i  -  T p2) 

=  (hrg-[g2  +  hCg-[g2)(Tg2  —  Tgi)  +  hcg2f(Tg2  —  Tf)  (214) 

For  the  packing  :  xeffapiS  =  hCptf(TPi  +  Tf) 

=  hrpig2(T  pi  +  rp2) 

+  drpib(Tpi  +  Tb)  (215) 

This  design  has  been  used  as  a  component  of  hybrid  heating 
system  [93,94].  To  simulate  the  performance  of  this  system,  a 
transient  analysis  has  been  proposed  by  the  authors  where  the 
above  energy  balance  has  been  written  for  a  rock  bed  solar  air 
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Fig.  17.  (a)  A  schematic  view  of  a  single  cover  solar  air  heater  (Choudhury  and  Garg  [88]).  (b)  A  schematic  view  of  a  single  cover  solar  air  heater  (Choudhury  et  al.  [88]). 


heater. 

dT 

MgiCgi  ^  —  otg\S-\-  hg-[g2(Tg2  —  Tgi)  —  hgia(Tgi  —  Ta)  (216) 
dT  i 

Mg2Cg2  ^  =  xg\ag2S  +  hpig2 (T pi  —  Tg 2)  —  hg2g\  (Tg 2  —  Tg  1) 

-  hg2a(Tg2  -  T f)  (217) 

dT 

M pi Cpi  —  Tgi  ^g2^pi*^  —  ^cpi f{T pi  ~~  T /) 

—  hrp-[g2  (T  pi  —  Tg2)  —  hpib(Tpi  —  T  jj)  (218) 

dT  f 

MfCf~dF=  ~mfcf(To  ~T fi)  +  ^p\f{Tp\  -  T f)  -hg2f(Tg2 

—  Tf)—  hcbf(Tb  -T f)  (219) 

Mf,C„ ^  =  hpW(T pl  -  Tb)  -  hcbf(Tb  -Tf)-Ub(Tb-  Ta)  (220) 


The  results  have  shown  that  the  hybrid  air-to-water  heating 
system  performs  much  better  when  coupled  with  a  rock  bed  air 
heater  than  when  coupled  to  a  conventional  empty  channel  air¬ 
heating  collector. 

2 A  6.  Double-pass  solar  air  heater  with  packed 

The  schematic  diagrams  of  the  double-pass  solar  air  heater  with 
packed  bed  are  shown  in  Fig.  18.  The  different  heat  transfer 
mechanisms  in  terms  of  the  various  heat  transfer  coefficients  are 
also  shown.  The  air  was  firstly  forced  through  the  packed  bed 
existing  in  the  upper  channel  (Fig.  18a),  formed  between  the  lower 
cover  and  the  absorber  plate,  and  is  then  re-circulated  to  flow  in 
the  opposite  direction  through  the  lower  channel,  formed  between 
the  absorber  and  back  plates.  Ramadan  et  al.  [96]  have  investigated 
heat  transfer  in  this  solar  heater.  The  mathematical  model  has  been 
developed  taking  account  of  the  following  assumptions:  (i)  the  air 
heater  operates  under  steady  state  conditions.  Flowever,  the  model 


remains  transient  due  to  the  time  dependence  of  solar  radiation 
intensity  and  ambient  temperature,  (ii)  The  heat  capacities  of  the 
glass  covers,  absorber  and  back  plates  and  insulation  were 
negligible,  (iii)  The  temperature  of  the  flowing  air  was  varied 
only  in  the  direction  of  flow,  (iv)  The  thermo  physical  properties  of 
the  flowing  air  were  assumed  to  vary  linearly  with  temperature 
(see  [7]).  (v)  There  was  no  temperature  gradient  across  the 
thickness  of  the  covers,  the  absorber  and  the  back  plates.  On  the 
basis  of  the  above  assumptions,  the  energy  balance  equations  for 
the  different  components  of  the  air  heater  were  written  as  follows: 


Cover 2  :  AgOtgS  —  Ag(hrgig2  +  flcgig2)(fg2  fgi) 

+  hwAg(Tg2  —  Ta)  +  Aghrg2Sby(Tg2  —  TSj<y ) 

(221) 

Cover  1  :  AgdgZgS  =  Ag(hrg-[g2  +  ficgig2)(7gi  —  ^g2) 

—  hrpigiAg(T  pi  —  Tgi)  —  AghCg-[  f\(T  —  Tg  i)  (222) 

Air  flow  in  upper  channel  :  Whp ^  (Tp  -Tfi) 

dT  n 

+  VWicpi/i  (Tpi  —  Tf- [)  =  rrif-[Cf  +  Whcg^f-i  (T^  —  Tg i) 

+  WUs(Tn-Ta)  (223) 

Absorber  plate  :  Ap/irplp(rpl  -  Tp)  =  Aphp/1(Tp  -  T fl) 

+  Aphpf2(T  p  —  T  f2)  +  Aphrpij(T  p  —  Tp)  (224) 

Air  flow  in  lower  channel  :  Whpf2(Tp  -T f2) 

+  wh f2b(Tb  -Tf2)  =  mf2Cf2-d^  +  WUs(Tf2-  Ta)  (225) 

Back  plate  :  Abhbf2(T f2  -  Tb)  +  Abhrpb(Tp  -  Tb) 

=  AbUb(Tb-Ta)  (226) 


The  various  heat  transfer  coefficients  given  in  the  above 
equations  were  calculated  using  the  correlation  given  in  the 
literature  [7].  A  suitable  computer  program  was  prepared,  for  the 
solution  of  the  energy  balance  equations  for  different  elements  of 
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Fig.  18.  (a)  Schematic  view  of  energy-flow  of  double-pass  type  solar  air  heater  with  recycle  (see  Ramadan  et  al.  [96]).  (b)  Schematic  view  of  energy-flow  of  double-pass  type 
solar  air  heater  with  recycle  (see  El-Sebaii  et  al.  [99]). 


the  solar  air  heater.  It  was  inferred  that  for  increasing  the  outlet 
temperature  of  the  flowing  air  after  sunset,  it  was  advisable  to  use 
the  packed  bed  materials  with  higher  masses  and  therefore  with 
low  porosities.  It  was  recommended  to  operate  the  system  with 
packed  bed  with  values  of  the  mass  flow  rate  equal  0.05  kg/s  or 
lower  to  have  a  lower  pressure  drop  across  the  system.  This 
mathematical  model  has  been  validated  by  comparisons  between 
experimental  and  theoretical  results  showing  that  good  agreement 
was  achieved. 

El-Sebaii  et  al.  [99]  have  proposed  an  investigation  of  heat 
transfer  on  a  double-pass  solar  air  heater  in  which  the  air  was 
firstly  forced  through  the  upper  channel,  formed  between  the 
lower  cover  and  the  absorber  plate,  and  then  re-circulated  to 
flow  in  the  opposite  direction  through  the  packed  bed  that  exists 
in  the  lower  channel,  formed  between  the  absorber  and  back 
plates  (Fig.  18b).  With  the  above  assumptions,  El-Sebaii  et  al. 
[99]  have  written  the  energy  balance  equations  for  the  system. 
For  cover  2,  energy  balance  equation  is  given  by  Eq.  (221).  For 
cover  1,  air  flowing  in  upper  channel,  absorber  plate,  air  flowing 
in  the  lower  channel  and  back  plate  may  be  written  as  follows, 
respectively: 

AgOtgXgS  —  Ag(hrg-ig2  +  frCgig2)(Tgi  —  Tg2) 

=  —  frrpigiAg(^pi  —  Tgi)  —  Aghcg-i  /i  (T/i  —  Tgi )  (227) 

Air  flow  in  upper  channel :  Whp^  (Tp  —  T^) 

dT  n 

+  Whcp\  f\(Tp\  -  I/O  =  mf-lCf-^-  +  WhcgA  /i(T/i  -  Tg i) 

+  WUs(Tfi-Ta)  (228) 


Absorber  plate  :  olvx\ApS  —  Aphcpf !  (T  p  —  T  ) 

+  Aphcpf2 (T p  —  T J2)  +  hp\iXpg\  (T p  —  Tgi )  +  Aphrpp\  (Tp  —  T p  1 ) 

(229) 

Air  flow  in  lower  channel  :  Whpf2{Tp  -T f2) 

+  Whcp ,  f2(T pi  -Tf2)  =  mf2Cf2-£-  +  WUs(T f2  -  Ta) 

+  Whbf2(Tf2-Tb)  (230) 

Back  plate  :  Abhbf2(T f2  -  Tb)  +  Abhrpib(Tpi  -  Tb) 

=  AbUb(Tb  -  Ta)  (231) 

The  various  heat  transfer  coefficients  given  in  the  above 
equations  were  calculated  using  the  correlation  given  in  literature 
[97,98].  The  program  developed  by  Ramadan  et  al.  [96]  was  used. 
On  the  basis  of  the  obtained  results  for  the  solar  air  heater,  it  had 
been  concluded  that  (i)  to  enhance  the  convective  heat  transfer 
coefficients  to  the  flowing  air,  as  well  as  the  heater  efficiency,  it 
was  advisable  to  operate  the  system  with  the  packed  bed  above  the 
absorber  plate  with  low  porosity  and  relatively  high  mass  flow 
rates  of  air.  (ii)  It  was  better  to  use  gravel  instead  of  limestone  as  a 
packed  bed  above  the  absorber  plate,  (iii)  It  was  recommended  to 
operate  the  system  with  and  without  the  packed  bed  with  values  of 
the  mass  flow  rate  equal  to  0.05  kg/s  or  lower  to  have  a  lower 
pressure  drop  across  the  system,  and  therefore,  reasonably  high 
thermo  hydraulic  efficiency,  (iv)  It  was  better  to  use  materials  with 
high  thermal  conductivity,  such  as  iron  scraps,  as  packed  beds 
above  the  absorber  plate  if  a  higher  outlet  temperature  of  air  is 
required  during  sunshine  hours.  Gravel  as  a  packed  bed  above  the 


1756 


R.  Tchinda /  Renewable  and  Sustainable  Energy  Reviews  13  (2009)  1734-1759 


Fig.  19.  Three  pass  two  covers  air  heater  with  in-built  thermal  storage  (Jain  and  Jain 
[100],  Jain  [101]). 


absorber  plate  during  low  intensity  solar  radiation  periods  was 
recommended. 

2.  I  7.  Multi-pass  solar  air  heater  with  in-built  thermal  storage 

Recently,  Jain  and  Jain  [100]  and  Jain  [101],  studied  an  indirect 
solar  crop  dryer  coupled  with  the  solar  air  heater,  investigated  the 
heat  transfer  in  the  multi-pass  solar  air  heater  with  in-build 
thermal  storage.  The  schematic  diagram  is  shown  in  Fig.  19.  In  this 
model,  the  solar  radiation  transmits  from  the  glass  covers  and  is 
absorbed  by  the  absorber  plate.  The  air  flows  in  between  the 
covers,  above  the  absorber  plate  and  below  the  storage  material, 
where  it  is  heated  along  the  path.  Jain  and  Jain  [100]  and  Jain  [101] 
have  written  the  energy  balance  equations  for  this  design  on  the 
base  of  the  following  assumptions  (i)  the  heat  capacities  of  the  air, 
covers,  absorber  plate  and  insulation  were  negligible;  (ii)  there  was 
no  temperature  gradient  along  the  thickness  of  cover;  (iii)  storage 
material  has  an  average  temperature  (Ts)  at  a  time  t,  (iv)  there  was 
no  stratification  existing  perpendicular  to  the  air  flow  in  ducts,  (v) 
the  system  was  perfectly  insulated  and  there  was  no  air  leakage. 
Based  on  these  assumptions,  energy  balance  equations  on  the 
double-cover  multi-pass  flat  plate  collector  with  in-built  thermal 
storage  were  given  as  follows: 

For  the  first  cover  :  ag-[Ag-[Sef  f  +  hrg2g\Ag\  (Jg2  ~Tg  i) 

=  ^cgl/l^gl  (Tgl  —  T/i)  +  frcglaAgl  (Tgl  — 

+  hrglskyAg l  (Tgi  —  Tsi(y)  (232) 

where  Tsky  =  Ta  -  6  (see  [102]): 

For  the  second  cover  :  otg2Tg\Ag2Sef f  +  hrpg2Ag2(Tp 
-  Tg 2)  +  hc f\g2Ag2 (T^i  -  Tg2)  =  hrg2g\Ag2 (Tg2  -  Tg  1) 

+  ^cg2  f2^g2(Jg2  -T  f2)  (233) 

For  the  absorber  plate  :  aprg  1  rg2Ag2Sef f  =  hrpg2Ap(Tp  -  Tg2 ) 


+  ^cp/2^p(^  p  —  T/2)  +  hcpsAp(T  p 

-Ts) 

(234) 

For  air  flow  in  the  first  channel  : 

Whcg  1  (Tgl  -  T  f-i ) 

=  mCf^L+whcfig2(Tn-Tg2) 

For  air  flow  in  the  second  channel 

:  Whcg2f2(Tg2  -  7 f2) 

(235) 

=  mCfd7^  +  Whcpf2(Tf,-Tg2) 

For  air  flow  in  the  third  channel  : 

Whcsf3(Ts-Tf3) 

(236) 

=  mCfd7ll  +  WUb(Tf3-Ta) 

(237) 

For  the  storage  material  :  hcps(Jp 

~  TS)AP 

=  MsCs^  +  hcsf3(Ts-Tf3)Ap 

(238) 

The  convective  heat  transfer  coefficients  from  the  plate  to  cover, 
parallel  to  each  other  and  inclined  at  the  angle  to  horizontal  has 
been  expressed  given  in  literature  [7,27]  and  in  the  appendix.  Jain 


and  Jain  [100]  have  computed  solar  intensity  on  the  inclined 
surface  by  using  the  method  given  by  Lui  and  Jordan  [103].  A 
computer  program  was  prepared  to  solve  the  above  set  of 
equations.  The  performance  evaluation  of  a  tilted  multi-pass  solar 
air  heater  with  in-built  thermal  storage  has  been  carried  out  for 
deep-bed  drying  applications.  The  thermal  energy  storage  also 
affect  during  the  off-sunshine  hours  is  very  pertinent  for  crop 
drying  applications.  With  the  above  mathematical  model,  Jain 
[101]  has  shown  that  it  is  useful  to  evaluate  the  thermal 
performance  of  a  flat  plate  solar  air  heater  for  the  crop  drying  in 
multiple  trays. 

2.18.  Matrix  (porous)  solar  air  heater 

The  porous  or  matrix  type  solar  air  heaters  are  shown  in  Fig.  20a 
and  b.  It  consists  of  single  glazing  at  the  top,  matrix,  back  metallic 
plate  and  insulation.  In  the  first  design  (Fig.  20a)  the  air  flow  is 
upward  through  the  matrix  and  in  the  second  design  (Fig.  20b),  it  is 
flow  downward  through  the  matrix.  These  models  have  been 
investigated  by  Sharma  et  al.  [98]  and  Sharma  et  al.  [104].  Sharma 
et  al.  [98]  analyzed  the  various  energy  exchanges  inside  the  porous 
bed.  The  energy  balance  equations  were  described  by  the  following 
equations  for  air  flowing  upward: 

f)T  f)2T 

MbCb  -^  =  kb-^-hfv(Tb-Tf)+  aaSexp(-ay)  (239) 

3T  r  d2Tf 

MfCf^f  =  GCf^yT-h^-Td  (240) 

The  corresponding  relations  for  air  flowing  downward  can  be 
obtained  from  the  above  equation  by  changing  the  sign  of  G.  The 
above  equations  were  solved  explicitly.  The  results  were  showed 
that  such  types  of  solar  air  heater  have  the  tendency  to  maintain  a 
constant  efficiency  throughout  the  operation.  The  experiment  was 
also  carried  out.  A  good  agreement  of  the  theory  with  experimental 
observations  almost  confirms  that  the  model  can  be  used  to  predict 
the  performance  of  the  porous  solar  air  heaters. 

A  transient  model  for  a  matrix  air  heater  having  a  single  glazing 
at  the  top  and  a  thin  metallic  plate  backed  by  an  insulator  at  the 
bottom  was  presented  by  Sharma  et  al.  [104].  The  mathematical 
model  used  was  similar  to  the  one  obtained  by  Sharma  et  al.  [98]  in 


Fig.  20.  (a)  Matrix  air  heater  (air  flow  upward)  (Sharma  et  al.  [98]).  (b)  Matrix  air 
heater  (air  flow  downward)  (Sharma  et  al.  [98,104]). 
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the  air  flow  downward  air  heater.  In  this  model  solar  insulation  and 
ambient  temperature  were  represented  by  Fourier  series  in  time 
with  y  dependent  Fourier  coefficients.  The  experimental  results  for 
a  varying  inlet  air  flow  rate  were  analyzed  and  the  rating 
parameters  were  included.  From  the  results,  it  was  observed  that 
the  matrix  air  heater  concept  shows  promise  for  collection 
enhancement  of  energy  at  a  cost  lower  than  most  of  other 
conventional  solar  air  heaters. 

3.  Conclusion 

Many  mathematical  approaches  have  been  developed  for 
simulating  the  performance  of  solar  air  heaters.  In  this  work,  a 
review  is  given  of  the  present  state  of  mathematical  modeling 
based  on  the  energy  analysis. 

The  models  on  energy  analysis  are  usually  expressed  in 
relatively  simple  equations  and  are  based  on  the  steady  state  of 


where  the  units  of  hw  and  V  are  W/m2  C  and  m/s,  respectively.  If  the 
unit  of  hw  is  changed  to  Kj/h  m2  K,  the  above  equation  becomes: 

hw  =  4.9  + 3.2  V  (A2) 

This  coefficient  has  also  been  given  by  Watmuff,  Charters  and 
Proctor  [109]  as: 

hw  =  2.8  +  3.0  V  (0  <  V  <  7  m/s)  (A3) 

An  empirical  equation  for  UT  was  developed  by  Klein  [67] 
following  the  basic  procedure  of  Hottel  and  Woertz  [68].  For  a 
horizontal  collector  with  two  glass  covers  UT  is  given  by  the 
relationship: 

UT  =  Un  +  Ur2  (A4) 

With 


[ _ 2(Tpm/520) _ I 

\[(Tpm/Tsky)/2  +  (1  +  0.089)h„  -  0.1166hw£p(l  +  0.7866  x  2)]0'43(1-100'tp»>)  J 

a(T  pm  +  TSky)(Tpm  +  7+) 

(fip  +  2  x  0.00591  hw)-1  +  [2  x  2  +  (1  +  0.089  hw  -  0.1166  hwsp)(l  +0.07866  x  2)  -  1  +  0.133  £p]/sp  -  2 


the  system.  Certain  assumptions  used  in  the  mathematical  form 
appear  to  be  too  ideal.  However  these  models  are  limited  by  the 
heat  transfer  efficiency  between  the  absorber  and  the  fluid.  Since 
heat  transfer  augmentation  usually  increases  friction,  optimum 
parameters,  such  as  mass  flow  and  collector  geometry,  have  to  be 
determined. 

Model  validation  is  a  key  step  in  model  development  since  it 
offers  the  possibility  of  comparing  computed  results  with  actual 
system  behavior.  Experiments  and/or  comparison  with  known 
analytical,  or  previous  numerical  results  are  mostly  used  to 
validate  the  mathematical  model.  Good  agreement  is  generally 
observed  and  the  computed  trends  are  also  compatible  with  those 
reported  for  experimental  observations.  However,  there  appears  to 
be  a  lack  of  agreement  among  authors  in  several  areas  and 
furthermore,  there  are  still  some  key  deficiencies  that  can  only  be 
addressed  through  continued  enhancement  of  the  existing 
modeling  techniques.  As  such,  there  exists  a  large  body  of  work 
to  predict  several  features  of  the  process. 

Future  work  should  be  directed  towards  a  more  comprehensive 
description  of  heat  transfer  phenomena,  which  would  involve  the 
generation  of  more  information  in  the  exergy  analysis  of  the  solar 
air  heaters. 
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Appendix  A 

The  convective  heat  transfer  coefficient  hw  for  air  flowing  over 
the  outside  surface  of  the  cover  depends  primarily  on  the  wind 
velocity  V.  MacAdams  [19]  obtained  the  experimental  result: 

(Al) 
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